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ABSTRACT 

Serendipitous observations of a pair z = 0.37 interacting galaxies (one hosting a quasar) show a 
massive gaseous bridge of material connecting the two objects. This bridge is photoionized by the 
quasar (QSO) revealing gas along the entire projected 38 ^pc sightline connecting the two galaxies. 
The emission lines that result give an unprecedented opportunity to study the merger process at this 
redshift. We determine the kinematics, ionization parameter (log C/ « —2.5 ± 0.03), column density 
{Nh,± « 10^1 cm-2), metallicity ([M/H]« -0.20 ±0.15), and mass (« IO^Mq) of the gaseous bridge. 
We simultaneously constrain properties of the QSO-host {Mdm > 8.8 x 10^^) and its companion 
galaxy {Mdm > 2.1 x lO^i; - 2 x IO^^Mq; stellar burst age=300 - 800 Myr; SFR~ 6 Mq 
yr~^; and metallicity 12 + log(0/H)=8.64 ± 0.2). The general properties of this system match the 
standard paradigm of a galaxy-galaxy merger caught between first and second passage while one of 
the galaxies hosts an active quasar. The companion galaxy lies in the so-called 'green valley', with a 
stellar population consistent with a recent starburst triggered during the first passage of the merger 
and has no detectable AGN activity. In addition to providing case-studies of quasars associated with 
galaxy mergers, quasar/galaxy pairs with QSO-photoionized tidal bridges such as this one offer unique 
insights into the galaxy properties while also distinguishing an important and inadequately understood 
phase of galaxy evolution. 

Subject headings: Quasars, Interacting Galaxies, AGN Feedback, Extended Emission Line Regions, 
QSO single: J204956.61-001201.7 



1. INTRODUCTION 

The paradigm of hierarchical structure formation is 
based on the frequent mergers of galaxies and their dark 
matter halos, as predicted and quantified by cosmolog- 
ical simulations which construct merger trees to track 
the merger history o f halos (e.g. iSpringel et al.l l2005bl : 
iFakhouri et al.|[2010( ). Obs ervations of galaxies 'caught 
in the act' of merging (e.g. lAT^ITOBfil : ILotz et al.l[2008b[ ) 
support this theory. Galaxy mergers are believed to be 
an important phase of galaxy evolution and are thought 
to be crucial to understanding super-massive black hole 
(SMBH) growth, galaxy morphologies, and the trunca- 
tion of star-formation within galaxies. 

While details of the merger process depend on 
the orbit, mass ratio and morphologies of the galax- 
ies involved, the basic sequence of events for a 

major (galaxy mas s ratio ^ 1:4) is rela tively 

well agreed upon (Mihos & Hernquist' '1 994 119961: 
Hopkins et al. 2008a; Cox ct al. 2006; Di M atteo et al.l 
20071: iChilingarian et al.ll2010( ). Galaxies under going a 



merger typically have orbital angular momentum that 
prevents direct collisions and leads to a series of close 
passages which have progressively smaller apogees due 
to dissipation in the merger. During their first close pas- 
sage, the galaxies exert strong tidal forces on each other. 
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For galaxies with a small bulge component, these tidal 
forces lead to bar formation. The stellar bar lags behind 
the gaseous bar inducing a torque on the gas which leads 
to the funneling of gas towards the center of the galaxy. 
The infiow also leads to higher central gas surface densi- 
ties which is predicted to spu r enhanced star formation 
(e.g. iMihos fc Hernquisti [19961 ). 

Subsequent close passages lead to smaller inflows as a 
large fraction of the available gas was used up during a 
first-passage 'starburst'. The second major phase of in- 
flow, therefore, occurs later due to a different mechanism 
as the galaxies are engaged in flnal coalescence. At that 
time, the remaining gas is brought rapidly to the cen- 
ter by dissipation into a dense compact central structure 
that rapidly forms stars. The final remnant galaxy in 
merger simulations is commonly an elliptical, dispersion- 
dominated system that has ceased star-formation be- 
cause its gas reservoir is used up in the starbursts and lost 
via stellar or AGN feedback (e.g. ISpringel et al][20 05a). 
However, the details of the effectiveness of 'quenching' 
star-formation and the actual process that is responsi- 
ble varies between different studies and depends on the 
choice of numerical prescription. 

For galaxies with large bulge components, the initial 
bar formation is suppressed and star formation is far 
less elevated in the early stages of the merger. Conse- 
quently, there is a large reservoir of gas remaining until 
the stages of final coalescence and the final starburst is 
more extreme. Because the galaxies coalesce over a short 
timescale, the final starburst can be more intense than 
that of a bulgeles s galaxy during first pa s sage. 

Some studies (jHopkins et al.l l2008at ISpringel et al.l 
I2005ai) argue that these merger events lead to elevated 
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levels of AGN activity and ultima tely a quasar (QSO) 
phase which is short (10^ — 10^-^ y rs; IMartini fc Weinberg! 
120011 : iMartini fc Schneided 120031 ) and extremely bright 
lO^^erg/s). This phase of the supermassive black 
hole (SMBH) evolution is thought to dominate the mass 
accretion history of the black hole (|Soltanlll982[ ) and can 
profoundly affect its host. During final coalescence, the 
large gas densities near the center of the galaxy (where 
the SMBH is located) leads to a large amount of efficient 
AGN fueling that manifests as a QSO phase. This sce- 
nario is supported by observations of quasar hosts show- 
ing disturbed mo rphologies (e.g. IBennert et al.l 120081 : 
iGreen et al.l |201C1[) consistent with being merger rem- 
nants. Observations disagree on whether mergers lead to 
elevated levels of AGN activity in the early stages of the 
merg er (i.e. between first an d second passage) with some 
(e.g., iWoods fc Gelled 120071) sug gesting an increase and 
others (e.g.. lEllisonet al. II2008D suggesting no increase 
in activity at these early stages. 

In addition to their obvious role in increasing galaxy 
mass, mergers provide a tantalizing solution to several 
observed cosmological trends. Firstly, they offer a con- 
venient explanation for the various SMBH-host rela- 
tions. Naively, the SMBH (which only has a gravitational 
sphere of influence « 1/1000 the size of the host) should 
have no correlation with the properties of its host. How- 
ever, SMBHs have been obser ved to be c orrelated with 
their host's st ellar mass (jM agorrian et al. 19 98), veloc- 
ity di spersion (jGebhardt et a l. 2000; Fcrrarcse fc Merritt] 
120001) , d ark matter halo mass (Ferraresei,2002i) , and mor- 
phology (jGraham et al.|[200TI ). If black holes are growing 
in galaxy mergers at the same time that many of these 
other galactic properties are changing, then this may ex- 
plain the observed correlations with the host. In fact, it 
has been proposed that this simultaneous evolution can 
lead to these variou s correlations even wit hout invoking 
any AGN feedback (ISutter fc Ricked [2010l) . 

Mergers may also be a key ingredient in quenching 
star-formation within galaxies. The basic picture is 
that galaxies can be broadly classified into one of two 
groups: (1) actively star forming spiral "blue cloud" 
galaxies or (2) passively evolving elliptical "dead" red se- 
quence galaxies. It appears that the red-sequence galax- 
ies s omehow evolve fr om actively star-forming galaxies 
(see iFaber et al.ll2007D . These hypotheses raise the ques- 
tion of how actively star-forming galaxies are "termi- 
nated" to produce red-sequence galaxies. While there 
hav e been many theories pr oposed (e.g., halo quench- 
ing; iDekel fc Birnboinill2006[ ). one of the more popular 
solutions is through the AGN activity, s tarbursts, and 
feedback induced by galaxy mergers (e.g. iHopkins et al.l 
l2008a ). In such a paradigm, the sequence of events (as 
described above) of the merger lead to a burst of star 
formation and AGN activity whose feedback can unbind 
a large fraction of the remaining gas. This quenches the 
galaxy while the merger also leaves the remnant with an 
elliptic al morpholog y, seemingly consistent with observa- 
tions (|Toomreill977( ). 

This theory is not without its flaws. Observational 
evidence suggests that galaxies transitioning from the 
"blue cloud" to red sequence (i.e. "green valley" galax- 
ies) do not appear to be dominated by merger remnants 
()Reichard et al.|[2009t ) nor the result of rapid quenching 



(|Martin et al.l [20071) . Nevertheless, it has established a 
theoretical framework that, in principle, can be tested 
with empirical observation. 

The details of SMBH fueling and feedback as well as 
star formation and supernova feedback are crucial to un- 
derstanding the actual role that mergers play in galaxy 
evolution. However, merger simulations cannot realisti- 
cally implement these processes due to the relevant phys- 
ical scales spanning a dynamic range of sub-pc to many 
kpc scales. Additionally, simulations are regularly un- 
able to accurately couple radiative transfer to the hydro- 
dynamics due to high computational expense. Simula- 
tors are thus forced to implement tuned sub-grid mod- 
els which sometimes are not physicall y well-motivated 
e.g., using Bondi-Hoyle accretion; see iBooth fc Schavi 
20091 ) ). Thus there remain many open questions about 
the processes at work during a galaxy merger that sim- 
ulations will be unable to enlighten for the foreseeable 
future. As such, observations are required to improve 
our understanding and guide future theoretical effort. 
However, observational studies of these processes are 
challenging for a number of reasons. In contrast to 
studies of nearby, merging galaxies which permit sen- 
sitive searches f or disturbe d morphologies and faint tidal 
structures fe.g.. IArp||1966[ ). such analysis at higher red- 
shift becomes prohibitive owing to the great distance. 
Low surface brightness tidal features are too faint to ob- 
serve and one requires much higher angular resolution 
to detect small morphological features at higher red- 
shift. Additionally, for all but the nearest galaxies (e.g., 
iHibbard et al.|[200ll ). observational studies of mergers are 
limited to studying the stellar distribution because the 
emission from the bulk of the gas is invisible at cosmo- 
logical distances. However, it has been shown that gas 
and stars may have completely different velocity struc- 
tures and distributions, especially in systems where a 
gala xy is undergoing a phase of enhanced AGN activ- 
ity ([Greene et al.l [20091) . Since the dynamics of the gas 
and not the stars ultimately fuel AGN activity and star 
formation, meaningful studies of these processes require 
the ability to study gas in mergers. Lastly, while the 
various feedback prescriptions may be tuned to match 
global properties of an ensemble population of galaxies, 
they must also be tested against actual merging systems. 

This gap in our understanding motivated a novel 
method for finding and studying merging galaxies that is 
inspired by the serendipitous discovery of a remarkable 
prototype. In our method, we search for systems where 
the ionizing flux of the quasar "lights up" the gas in a 
merger in order to find new observational constraints on 
merger models. 

With a n ionizing photo n flux (S) of ~ 10^^ pho- 
tons/s/sr dTadhuntedll996[ ). a quasar can ionize gas to 
conditions comparable to H II regions out to distances 
of 60 kpc. While observations of highly ionized gas 
with large (tens of kpc) spatial extent known as ex- 
tended emis sion line regions (EELRs) have been known 
for decades (Matthews fc Sandagc 1963; Wamolcr et aU 
19751: fStockton fcJMacKentv 1987; Penston ct al. 199^^ 
Fu fc StocktonI [20091 ). our study focuses on a novel ap- 
plication of this standard result where we utilize quasar 
photoionization to study a pair of merging galaxies. In 
detail, we focus on an EELR produced by a quasar that 
has been triggered by a merger event. During such a 
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merger, one also expects extended tidal features of gas 
out to distances of 100 kpc that may be photoionized by 
the quasar and could thus be rendered visible through 
recombination and forbidden lines. Observations have 
confirmed that such extended material that may be as- 
sociated with a galaxy merg er can be seen in emissio n 
lines near at least one quasar (|Villar-Martm et al.ll2Q10l ). 

Our prototype was serendipitously discovered while se- 
lecting targets for a study of gas in the extended ha- 
los of galaxies via quasar sight lines (Tumlinson et al 
in prep.). That study used the Sloan Digital Sky Survey 
DR5 (iAdelman-McCarthv et al.ir2067l) to find foreground 
galaxies close to quasar sight lines based on photomet- 
ric redshifts. However the galaxy they had hoped to be 
in front of J2049-0012 (a z=0.369 quasar) by a cosmo- 
logical distance was in fact found to have the same red- 
shift. Further inspection revealed that the long-slit spec- 
trum of the system contained bright emission lines along 
a remarkable gas "bridge" connecting the quasai0 to its 
companion galaxy, 38 kpc away in projection, and 
photoionized by the hard spectrum of the quasar. This 
bridge is most likely the result of a tidal interaction of the 
galaxies during the first passage of their ongoing merger. 
The extended gas around these objects, normally com- 
pletely invisible at this distance, is being lit-up by ra- 
diation from the quasar. This provides an unparalleled 
opportunity to study the kinematics and structure of the 
gas in the merger, which otherwise we would not have 
even known was occuring. At the same time, we have 
detailed measurements of the galaxy properties, includ- 
ing the outer regions of the quasar's host galaxy. The 
companion galaxy is also a member of the 'green valley' 
and we are afforded the chance to examine this poorly 
understood and rare phase of galaxy evolution in this 
merger context. 

This is the first in a pair of papers on this particular 
system. This paper focuses on the details of the merging 
system (mass ratio, timescales, separation, merger stage, 
companion galaxy properties, SMBH mass and luminos- 
ity) as well as discussion of the emission line bridge and 
an explanation for its ionization source. Paper II focuses 
on the quasar including its lifetime, isotropy, impact on 
companion galaxy, implications for its triggering, and the 
placement of the SMBH on Magorrian-like relations at a 
peculiar phase of its evolution. 

We assume a cosmology with f^A — 0.73, Hm = 0.27, 
and Ho = 72 /172km s"^ Mpc"^ 

The paper is organized as follows: § [2] describes the 
observations of this system; § [3] gives an overview of the 
results of the observations; § 4 details our measurements 
of emission line fluxes and kinematics; § 5 describes the 
inferred properties of the interacting pair of galaxies; § 6 
discusses our analysis of the bridge connecting the galax- 
ies including our analysis of its origin; § 7 presents further 
discussion; § 8 presents a summary. 

2. OBSERVATIONS & DATA REDUCTION 

^ We reserve the terms 'quasar' and 'host galaxy' to refer to 
the galaxy that is undergoing a quasar phase while terms with the 
adjective 'companion' refer to the other galaxy with which the QSO 
host is interacting. We interchange the use of the terms QSO and 
quasar and make no distinction between the two classes of objects 
for the purposes of this paper. 



We obtained a set of spectroscopic and imaging ob- 
servations of this quasar /galaxy pair using the dual- 
channel Low R esolution Imaging Spectrometer (LRIS; 
lOke et aP 119951 ) on the Keck I 1 0m telescope and the 
dual-channel Kast spectrometer ([Miller fc Stond Il993f ) 
on the 3m Shane telescope at Lick Observatory (see Ta- 
ble H]). 

Two long slit spectra were taken using LRIS with a 1.0" 
slit width, D560 dichroic, 600/4000 grism, and 600/7500 
grating tilted to an angle of 28.15°. This gave a red 
wavelength coverage of roughly 5600 to 8200A, a dis- 
persion of 0.63 A/pixel, and FWHM « 160 km s^^. 
In the blue we had coverage of 3100 to 5600 A with 
FWHM ss 240 km s~^. Subsequently we took further 
LRIS observations with the grating angle tilted to 29.15 
giving us coverage of 6300-9500 A to cover Ha and [N II] 
emission of our z=0.37 system. We also took a higher 
resolution spectrum covering the [O III] emission lines 
to better characterize the gas kinematics. We used the 
0.7" long slit, D560 dichroic, and 1200/7500 grating 
blazed to an angle of 40.15. This gave a dispersion of 
0.4 A/pix~^, FWHM « 65 km s~^, and wavelength cov- 
erage of roughly 6100 - 7000A. The FWHM for each of 
these setups was estimated by measuring the parameters 
of Gaussian models fit to sky lines in the science expo- 
sures. 

All of the above spectra were taken with the long slit 
aligned along the line connecting the two sources (so that 
both the quasar and companion galaxy were within the 
slit). Another spectrum was taken using the Lick Ob- 
servatory's Shane 3-meter with the Kast spectrograph. 
This spectrum was taken with a position angle (PA) per- 
pendicular to the line connecting the quasar/galaxy pair, 
roughly halfway in between (such that neither the quasar 
nor the companion were in the slit). This spectrum was 
taken with the d55 dichroic and 600/7500 grating. This 
grating was tilted to span a wavelength range of approx- 
imately 5500 - 8300A with a dispersion of 2.37 A/pixel. 

Additionally we imaged the system with the LRIS cam- 
era for one 60, one 200, and two 230 second exposures 
in the B band as well as one 90, and three 200 second 
exposures in the R band (2008-10-03). The seeing was 
« 0.64" FHWM. 

We reduced the spectra using the Low-Redu}(0 
pipeline developed by J. Hennawi, S. Buries, D. Schlegel, 
and J. X. Prochaska. We implemented this program to 
perform the following procedures: (1) process the flats, 

(2) use arcs to make a 2-dimensional wavelength image, 

(3) make a slit profile, (4) process the images, (5) identify 
objects in the slit and trace them, (6) perform sky sub- 
traction, (7) extract the spectra, (8) correct for flexure of 
the spectrograph using sky lines, and (9) create a sensi- 
tivity function from a standard star. Slight modifications 
to the standard pipeline were necessary to perform sky 
subtraction due to the extended emission and the prox- 
imity of the galaxy to the quasar. Specifically, we turned 
off local sky subtraction after verifying that a global fit 
provided a good estimate of the sky background in the 
region of interest. Next, we coadded the images using 
Long-Coadd2d, a code developed by Robert da Silva 
and Michele Fumagalli, removing cosmic rays through 

^ http:/ /www. ucolick.org/~xavier/LowRedux/ 
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TABLE 1 

Observations 



Instrument 


Slit 


Dichroic 


Camera 


Disperser'^ 


Resolution [km/s] 


Filter 


Exposure 


Standard Star 


Date 




["] 










(UT) 




Keck/LRIS 


1.0 


D560 


R 


600/7500 (28.15) 


160 




900x2 


Feige 110=" 


2008-10-03 


Keck/LRIS 


1.0 


D560 


B 


600/4000 


240 




830x2 


Feigc 110=^ 


2008-10-03 


Keck/LRIS 


1.0 


D560 


R 


600/7500 (29.15) 


160 




300x3 


BD+28'' 


2009-09-17 


Keck/LRIS 


0.7 


D560 


R 


1200/7500 (40.15) 


65 




600x2 


BD+28'' 


2009-09-17 


Shane/Kast 


1.0 


d55 


R 


600/7500 (10775.6) 






1800 




2009-05-26 


Keck/LRIS 




D560 


B 






B 


60,200,230x2 




2008-10-03 


Keck/LRIS 




D560 


R 






R 


90,200x3 




2008-10-03 



^ (Bohlin 1996) 

^ { Bohli n ct al. 200ll 'l 

Values in parenthesis denote the grating angle. 

Measured as the FWHM of a sky line. 

comparisons of the separate exposures. 

The atmosphere and instrument response influence the 
intensity and shape of the recorded spectrum. For point 
sources, the standard procedure to convert the observed 
spectrum into astronomical flux units is to construct a 
sensitivity function from the spectrum of a spectrophoto- 
metric standard star observed with the same instrumen- 
tal configuration. This is done in 1-dimension, meaning 
that a 1-dimensional spectrum is extracted from the 2D 
output spectrum from the instrument and is then com- 
pared with a ID model. We performed this procedure 
and applied the resulting sensitivity function to the ob- 
jects extracted in the slit. 

Unfortunately the region near the dichroic (A « 5400 — 
5600A) was particularly difficult to calibrate for the 2008- 
10-03 observations. In particular, we find that the flux- 
ing resulted in a galaxy spectrum that had an unphysical 
behavior near the dichroic. We correct for this through 
comparison of our Keck/LRIS spectrum of the QSO and 
the SDSS spectrum of the same QSO. Although the ab- 
solute and relative fluxes of quasars are known to vary, 
the color variation is generally small for rest wavelengths 
greater than 2500 A ( Wilhite et al.ll2(305iV Avoiding is- 
sues that may arise from Fe emission lines we only sought 
to flt the general shape of the continuum across the 
boundary. In order to account for any discrepancies be- 
tween ours and the SDSS fluxing, we matched the spectra 
by allowing our Keck/LRIS spectrum to be multiplied by 
a polynomial in wavelength. To find the best such 2nd 
order polynomial we used the singular valued decompo- 
sition method. Once the overall shapes of the spectra 
was matched using the above procedure, we fit a 3rd 
degree polynomial to each side of the dichroic for both 
the SDSS and matched LRIS 2008-10-03 spectra. The 
ratio of those polynomial fits was used to correct the 
LRIS spectra. Thus the spectrum correction forces the 
LRIS spectrum of the QSO to have the same shape in 
the region around the dichroic as the SDSS spectrum. 
This same correction was then applied to the companion 
galaxy spectrum. 

The final product from the above procedure was a sen- 
sitivity function that converts electrons per second per 
spectral pixel into astronomical flux units. This allows 
us to calibrate the extracted, ID spectrum of any object 
in the slit. Our data, however, features extended emis- 
sion that spans many arcseconds spatially. Therefore we 



also require a 2-dimensional flux calibration that con- 
verts electrons per second per detector pixel into surface 
brightness units. To accomplish this task, we needed to 
apply the 1-dimensional fluxing (which has a mapping 
of pixel to sensitivity that is accurate at the center of 
the slit), to the entire 2D spectrum. Assuming that the 
appropriate correction factor is only a function of wave- 
length (specifically that there was a 1-to-l mapping of 
wavelength to sensitivity) we can apply this correction 
to the entire CCD using the 2-dimensional wavelength 
image constructed with Low-Redux from the arc-lamp 
images. Interpolation of our sensitivity function evalu- 
ated at every pixel's wavelength gives a correction factor 
which we then apply to the spectra. 

Since we had taken long slit data, each of our spectra 
had portions of the CCD illuminated only by the sky. We 
used these to test the 1-to-l assumption described above 
by measuring sigma-clippcd medians along the spectral 
direction across spatial regions. The spatial variation in 
sensitivity along the slit was found to be at most a 1% 
effect. 

3. J2049-0012: OBSERVATIONAL CHARACTERISTICS 

This section provides a qualitative overview of the 
properties of the J2049-0012 system. More quantitative 
analysis is presented in the sections that follow. 

3.1. Imaging 

This quasar/gal axy pair was selected from the SDSS 
Data Release 5 (lAdelman-McCarthv et al.l |2007| ) and 
GALEX DR3 (Martin et al.' '2005D catalogs as a candi- 
date system for an HST/COS survey designed to study 
the gas surrounding z ^ 0.2 galaxies with absorption- 
line spectra of UV bright, background quasars (Cy- 
cle 17, ID=11598; PI: Tumlinson). The quasar pho- 
tometry and redshift {zq = 0.3693) were known from 
SDSS observations while the galaxy's SDSS photome- 
t ry yields a photom etric redshift Zphot = 0.28 ± 0.09 
(|Qvaizu et alj 120081) . suggesting it lay foreground to 
the quasar. Fig. [1] presents the R band imaging from 
Keck/LRIS. The quasar (J2049-0012) is the bright 
point source at the center of the image and the galaxy 
(J204956.61-001201.7) is offset by 8" (38 h:f^ projected 
kpc) to the SW (lower right). 

The companion galaxy as seen in an SDSS find- 



5 



TABLE 2 

Compiled Photometric Measurements 



Quasar 


Galaxy 


Extinction^ Ref 


(mag) 


(mag) 


(mag) 



FUV 


18.54 ±0.108 




0.46 


1 


NUV 


18.44 ±0.063 




0.49 


1 


u 


18.05 ±0.01 


21.42 ± 0.28 


0.46 


2 


9 


17.89 ±0.01 


20.93 ± 0.05 


0.34 


2 


r 


17.79 ±0.01 


19.87 ±0.03 


0.25 


2 


i 


17.75 ±0.01 


19.62 ± 0.04 


0.19 


2 


z 


17.10 ±0.01 


19.46 ± 0.11 


0.13 


2 


J 


16.51 ±0.142 




0.08 


3 


H 


15.74 ±0.150 




0.05 


3 


K 


14.62 ±0.101 




0.03 


3 



Calculated using Ay = 0.292 from ISchlegel et al. l l fT99l) 
and a ssuming the dust extinction law from [Cardclli c t al] 
II1989I V All magnitudes listed have not been extinction cor- 

r ected. 

1 ITrammell et al.| 120071 ') 
2 [Xbazaiian ct al. ( 20091) 
3 [UTTtri et al., (i2003i) 





Fig. 1. — R-band image (which covers the redshifted [O III] and 
H/3 lines) zoomed in on the position of the galaxy /quasar pair. 
The image is 17" on a side and is oriented so that N is up and E 
is left. One can discern the bridge connecting the two galaxies as 
well as some possible tidal debris. The quasar (Q) and companion 



galaxy (G) are separated by 8" 
z = 0.3693. 



corresponding to Z&hj,, kpc at 



ing charlQ is a peculiar green, uncommon for galax- 
ies observed in the SDSS. This is illustrated in Fig. [2] 
which compares the restframe color and absolute mag- 
nitude of this gala xy to a sample of DEEP2 galaxies 
at 3 < z < 0.4 foavis et al.l I200I [2007L iCoil et alJ 
12001 IWillmer et al.ll200fil ). We find that the color of 
J204956. 61— 001201.7 places it between the two primary 
populations of red and blue galaxies. Such galaxies are 
commonly described as lying within the 'green valley', a 
region that possibly marks the transition fro m blue star- 
formi ng galaxies to red, quiescent systems (jFaber et al.l 
[20071) . The processes that drive galaxies into the green 



http: / /cas.sdss.org/dr7/en/tools/chart/chart.asp?ra= 
0.20136 



GQ 




Fig. 2. — A rest-frame color-magnitude diagram of galaxies in 
the redshift slice 0.3 < z < 0.4 from the DEEP2 survey (grey 
points; [Coil ct al. 2004). The companion galaxy of the J2049-0012 
system (z=0.37) is represented by a black diamond. It appears in 
the 'green valley" between the main loci of the blue cloud and red 
sequences. The dotted line marks the sep aration between the r ed 
sequence and blue cloud as established bv lWillmer et aT] 120061) . 



valley are not well established. Scenarios include merger- 
induced star bursts where the tidal forces the galax- 
ies e xert on each other f unnel gas towards their cen- 
ters (|Quintero et al.l [20041) . which in turn fuels a mas- 
sive starburst and u ltimately a bright quasar phase 
([Hopkins et al.' 2008b). Such a model predicts that the 
green valley region should coincide with AGN activity, 
which is supported by recent observational e vidence (e.g. 
iSchawinski et al.l I20Q9I: iNandra et all I2007D . However, 
the level of AGN activity in the companion galaxy in 
our quasar/galaxy pair is unclear. The Id extraction 
of our galaxy shows faint [Ne V] AA3426, 3346 emission 
(a commonly used indicator of AGN activity) , but even 
this emission may be the result of photoionization by the 
neighboring quasar. 

Additionally, such a model may also predict that ob- 
jects in the 'green valley' have transitional stellar popu- 
lations such as those of post-starburst ga laxies. This is 
supp o rted by observationa l studies (e.g., iVergani et al.l 
l2010t iKocevski et al.l I2009D . These galaxies have their 
light dominated by a stellar population that has aged 
past the point where its stellar light is dominated by O 
and B stars. Examination of our spectra (^ 13.2^ . suggests 
that our galaxy may have such a stellar population with 
some ongoing star formation. 

An alternative explanation for the companion galaxy's 
location in the 'green valley' may be that it is a blue star 
forming galaxy seen nearly edge on such that dust makes 
it app ear redder than it would otherwise (Brammer et sA\ 
120091 ). As noted in the following section, there is evidence 
of strong rotation in the galaxy which may suggest that 
the galaxy is inclined, however there is little evidence for 
strong dust absorption from the observed Ha/H/3 ratio. 

Morphological analysis of the two objects is difficult. 
Firstly, the quasar (combined with its host galaxy) is 
consistent with a point source. The companion galaxy, 
meanwhile, is resolved with only a few spatial elements 
(observed to have FWHM of 1.86" [9.5 kpc] when our 
seeing was w 0.64") and is als o consistent with h aving 
3l2.48392&decaai axis ratio of 1. Following iLotz et all ([2008aD and 
iConsehca (|20()3[) . we measure a concentration for the 
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Fig. 3. — One dimensional Keck/LRIS spectra of the galaxy and quasar comprising our interacting pair (the J2049— 0012 system). Note 
the strong Balmer series absorption lines in the galaxy spectrum; these are suggestive of a starbust population. Also note that the quasar 
and galaxy are at nearly identical redshift. The quasar is at = 0.3691 ± 0.001 and the companion galaxy is offset by only 160 ± 20 
km s~^. As a clarification, we note that this quasar velocity is not the systemic velocity of its host galaxy. We find the veloc ity difference 
between the two galaxies to be 30 ± 30 km/s. The strongest absorption and emisson lines are labelled by their ion. Fig. 1101 shows an 
enlarged view of the post-starburst features along with models used to constrain the age of the stellar population. 

companion galaxy of C = Slogj^g {^~^ where rso is the 
radius containing 80 percent of the total light, r2o is sim- 



ilarly the radius containing 20 percent of the total light, 
and we have defined the total light to be the light within 
1.5 times the Petrosian radius. We find that the galaxy 
has a concentration of 2.61 classifying it as a spiral (the 
standard dividing line is spirals have C < 3). However, 
we note that since viq is within one seeing element, our 
measured concentration is only a lower limit and there- 
fore the profile may actually be more consistent with an 
elliptical or bulge dominated system with a higher con- 
centration. Around both galaxies there is some hint in 
the deeper R-band exposure of tidal debris, but since 
those regions did not fall into our slit we have no spec- 
troscopic information about them and it is difficult to 
determine if this emission is truly associated with the 
galaxy or only close in projection (see Fig. [Ij. 

Table [2] lists photometry for the two objects compiled 
from a variety of surveys, spanning UV to near-IR fre- 
quencies. 

Examination of the photometric redshifts of nearby 
galaxies reveals no obvious overdensity at the quasar's 
redshift. We examined this by looking at the number of 
objects within 1 Mpc projected distance from the QSO 
with consistent photometric redshifts (1.5ct) measured 
for 1000 random positions in the SDSS footprint that 



contained at least a single galaxy (at any redshift). The 
quasar appears in the low end of this distribution and 
hence there is no evidence for the quasar being in a rich 
and large cluster. 



3.2. ID Spectra of the QSO Host and Companion 

The discovery spectrum was taken with the origi- 
nal goal of confirming this system as a projected back- 
ground/foreground, quasar/galaxy system separated by 
a cosmological distance. We observed the objects us- 
ing the Keck/LRIS spectrometer with its 1" long slit 
aligned across the pair. Fig. [3] shows the fluxed ID spec- 
tra. The quasar shows strong and broad H/3 and Mg II 
emission lines confirming its SDSS classification and red- 
shift. A Gaussian fit to the [O III] emission lines gives 
Zq = 0.3691 ± 0.001 which is within the errors of the 
SDSS measured redshift of 0.3693 ± 0.0015. The broad 
emission lines appear unusually asymmetric. 

The companion galaxy's photometric redshift placed it 
well in front of the quasar, however the spectrum reveals 
the galaxy has strong emission and absorption features 
at nearly the identical redshift of the quasar with only 
a (5u ss 150 km/s second offse10 (Fig. [3]). The spectrum 

^ This velocity difference, however, does not precisely reflect the 
veloc ity difference between the two galaxies which is discussed in 
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Fig. 4. — 2D discovery spectrum of the J2049-0012 quasar/galaxy system centered near H/3 and [O III] emission made using the 
Keck/LRIS 600/7500 grating. The continuum light of the quasar (top) and galaxy (bottom) are the bright horizontal bands in the image. 
The bright quasar spectrum is centered at the kpc (O") spatial position and saturates the image in this stretch. The image covers, from 
left to right, the H/3, [O III]A4959, and [O III]A5007 transitions corresponding to 2 = 0.37. These emission lines are detected for the quasar, 
its host galaxy, and the companion galaxy as well as a 'bridge' of emission that connects the two objects. The emission lines are not 
detected beyond the spatial slice presented here. 

3.3. 2D Spectrum 
In Fig.m we present a slice of the co- added, 2D discov- 



is notable for strong emission lines characteristic of pho- 
toionization together with strong Balmer series absorp- 
tion lines that indicate a moderate-aged stellar popula- 
tion (tagc ^ 100 Myr). The absorption lines may suggest 
that the galaxy has recently undergone a 'burst' of star- 
formation that may now be fading. After O and B stars 
die after ~100 Myr, A stars remain a significant contrib- 
utor to the stellar light until their death approximately 
1 Gyr later. Type A stars provide characterisitic Balmer 
series absorption lines that make them easy to distin- 
guish from other stellar types. If there is an older stellar 
population (t > 1 Gyr) contributing to the galaxy light, 
the K giants will show strong Ca H-l-K absorption lines 
(jPressler fc Gunnlll983l) . The galactic spectral type that 
is characterized by both Balmer absorption lines and Ca 
H-fK lines and with no emission lines is commonly re- 
ferred to as K-t-A post-starburst, a re latively rare phas e 
for galaxies especially at low redshift (jWild et al.ll2009( ). 
These indicators suggest that the galaxy's star formation 
rate is fading and/or has recently shut off. A common 
scenario used to explain such a signature is a passively 
evolving stellar population that underw ent a burst of star 
form ation that has recently subsided (jPressler fc Gunnl 
I1983D . Such bursts may be caused by tidal interactions 
in a merging galaxy sys tem (Snyder et al. in prep.; 
iMihos fc Herngulstl 119960 . This applicability of this in- 
terpretation for our system is further discussed in i)5.2.3l 
Because the galaxies are at nearly the same redshift, 
there is little information that can be used to determine 
which galaxy is foreground. Mg II absorption is clearly 
visible in the companion galaxy's spectrum, while none is 
evident in the spectrum of the quasar (See Fig. ^ . This 
could indicate that the companion galaxy is background 
to the quasar host galaxy (e.g. Rubin et al. 2010). How- 
ever, the similar velocities of the two galaxies precludes 
one from unambiguously associating the gas with either 
the quasar host or the companion galaxy. 



ery spectrum of the system centered near the observed 
II/3 and [O III] emission lines. There are two sources with 
visible continua which appear as horizontal lines in this 
2-dimensional spectrum: the quasar and the companion 
galaxy. The [O III] and H/3 lines exhibit faint emission 
oriented at « 20° counter clockwise from vertical that ex- 
tends « 2" (10 h^^ kpc) away from the spatial location 
of the quasar in directions both towards and away from 
the companion galaxy. We associate this emission with 
the underlying host galaxy of the quasar. The weaker 
continuum source located at impact parameter p w 8" 
(38 kpc) and its associated emission lines are from 
the companion galaxy. The galaxy emission shows obvi- 
ous signatures of rotation; the [O III] emission spans ~5 
A corresponding to Av — cAX/X « 300 km s~^ which 
implies a rotation speed Vrot > 100km s^^. Interestingly, 
the galaxy's emission properties vary as a function of spa- 
tial distance from the quasar: H/3 emission is stronger on 
the side away from the quasar than the side toward it, 
while [O III] shows the opposite trend. 

Of greatest interest to this work is the detection of a 
'bridge' of emission that connects the quasar and galaxy. 
The bridge is strongest in [O III] but H/3 and other 
lines shown in Fig. |6] are also positively detected. It 
appears that the bridge lies roughly along a line con- 
necting the spatial and velocity centers of the two galax- 
ies. Fig. H] also reveals that along the line connecting 
the quasar and galaxy, the emission is confined to the 
space between them. In this respect, the emission dif- 
fers from the EELRs comm only observed around quasars 
(e.g. iFu fc Stocktonll2009f ). i.e., it appears to be directly 
related with the gas in this interacting system. To ex- 
plore this point further, we obtained an optical spectrum 
with the Kast spectrometer on the 3m Shane telescope 




I 

Fig. 5. — A cartoon depiction for the assumed geometry of the 
J2049— 0012 system. We assume that the bridge forms a cyhnder 
of diameter D connecting the quasar host galaxy to the compan- 
ion galaxy, ri and are radial distances from the quasar. The 
volume density of hydrogen is assumed to be a constant nj^ . Also 
noted are the approximate slit positions of both the Keck/LRIS 
and Shane/Kast spectra. Using this simple model, we can estimate 
properties of the bridge including its volume and column densities 
and its mass. 

at Lick Observatory with the long sht centered on the 
bridge but oriented perpendicular at a PA perpendicular 
to that of our Keck/LRIS spectrum. By fitting a gaus- 
sian to the Kast spectrum, we found its FWHM=2.0", 
which was consistent with our seeing limit for the night. 
Therefore we constrain the spatial extent of the bridge to 
be < 10.2 K^^ kpc perpendicular to the line connecting 
the quasar and companion galaxy. 

3.4. Summary 

Fig. [5] presents a cartoon of the inferred layout of the 
system. Taken altogether, the observations of this sys- 
tem suggest the following scenario: we are observing 
a pair of merging galaxies that have undergone an in- 
teraction which has tidally stripped material from at 
least one of the galaxies. This is suggested by the 
nearly identical redshifts of the galaxies and the pres- 
ence of the bridge connecting the two galaxies. Such 
material is a standard predic tion in theoretical models 
of gas-rich galaxy encoun ters (jToomre fc Toomrill972al : 
IMihos fc Hernguistl I1996D as well as a commonly ob- 
served phenomenon in nearby, interacting galaxies (e.g. 
lA^ [19661 iHibbard et all mm . The ongoing merger 
could have triggered a burst of star-formation in one (or 
both) of the galaxies and AGN activity in the galaxy 
now observed to be in a quasar phase. The companion 
galaxy appears to be in the 'green valley' and has a possi- 
bly fading starburst stellar population. This population 
is consistent with a starburst that was caused during the 
first passage of a galaxy interaction and may indicate 
that the galaxy is likely undergoing a substantial change 
in its stellar light properties. The fact that the higher 
ionization state lines are brighter closer to the quasar 
(i i4.2[) suggests that the quasar is shining on the tidal 
material, photoionizing the majority of that gas. If this 
hypothesis is correct, then the full toolbox of low den- 
sity astrophysics allows us to interpret observations of 
the bridge's emission lines to give insights into the tidal 
material's dynamics, column density, mass, temperature, 
metallicity, and volume density. In the following sections 
we test this physical interpretation of the system and its 
implications. 

4. EMISSION LINE MEASUREMENTS 
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Fig. 6. — Portions of the discovery spectrum taken with 
Keck/LRIS on 2008-10-03 after continuum subtraction. The star 
denotes the spatial position of the quasar and the open circle marks 
the position of the companion galaxy. A grey bar has been placed 
over the core of the quasar to aid the eye. These images have 
been interpolated to be on the same spatial and velocity grid. The 
[O III]A4959 panel shows the spatial boxes used throughout the pa- 
per which separates regions in each of the bridge (Ob, lb, 2b) and 
the companion galaxy (3g,4g). 

The flux measurements and line ratios constrain prop- 
erties of the gas and galaxies involved in this system. 
The kinematics provide estimates for the masses of the 
galaxies and give information about the merger stage and 
geometry. We present the emission line fluxes in Table |3l 
Inspection of Fig. |6] reveals that the emission line fluxes 
and the line ratios vary along the slit hence we have bro- 
ken the bridge into three bins (Ob, lb, and 2b) and the 
galaxy into two equally sized bins (3g and 4g) as labelled 
in Fig. [HI The bounds of the regions noted as "bridge" 
and "galaxy" are determined by eye through inspection 
of the kinematic profiles. In the higher resolution spec- 
trum, it is clear that the bridge has distinct kinematic 
characteristics from either galaxy. 

4.1. Line Fluxes 

To properly isolate the emission line fluxes, we needed 
to flrst subtract the continuum light from the quasar 
and the galaxy. This required construction of an ac- 
curate model of the flux from each object. We took the 
2-dimensional proflle flts from Low- Redux (interpolat- 
ing over the emission-line regions) and multiplied by the 
1-dimensional boxcar extraction. An additional compli- 
cation arises from the companion galaxy's Balmer ab- 
sorption lines. To account for these spectral features, we 
used the stellar population model described in § 15.2.31 
to model the absorption for the Ha and II/3 fluxes. This 
resulted in a spectrum containing only the emission lines 
(Fig.©. 

The main contributions to the error are the statistical 
errors of the signal (Poisson errors) and the errors asso- 




TABLE 3 

Emission Line Measurements 





Region^ 


Flux 
10^-^^ erg/s/cm^ 


Dereddened Flux^ 
10"-'^'' erg/s/cm^ 


Flux c 


[N II]A6585 


Ob 
lb 
2b 

3g 

4g 


<0.87d 

I. 34±0.38 
<1.29 

II. 68±0.60 
13.42±0.67 


<0.99 
1.52±0.43 

<1.46 
13.28±0.68 
15.26±0.76 


<0.07 
0.12±0.03 

<0.05 
0.27±0.01 
0.75±0.04 


Ho 


Ob 
lb 
2b 
3g 

4g 


4.89±0.43 
3.53±0.50 
5.19±0.56 
32.51±0.77 
38.18±0.79 


5.56±0.49 
4.02±0.57 
5.90±0.64 
36.98±0.88 
43.44±0.89 


0.38±0.04 
0.31±0.04 
0.22±0.02 
0.77it0.02 
2.14±0.08 


[N II]A6549 


Ob 
lb 
2b 
3g 
4g 


<0.92 
<1.18 
<1.32 
4.48±0.60 
4.95±0.64 


<1.04 
<1.35 
<1.50 
5.10±0.68 
5.64±0.73 


<0.07 
<0.10 
<0.06 
O.llitO.Ol 
0.28±0.04 


[O III]A5007 


Ob 
lb 
2b 
3g 
4g 


11.87±0.30 
10.65±0.35 
21.85±0.41 
39.24±0.55 
16.52±0.53 


14.62±0.37 
13.11±0.44 
26.90±0.51 
48.30±0.67 
20.34±0.65 




[0 III]A4959 


Ob 
lb 
2b 

3g 
4g 


4.27±0.26 
3.70±0.31 
7.41±0.36 
12.93±0.48 
4.77±0.46 


5.27±0.32 
4.57±0.39 
9.14±0.44 
15.96±0.60 
5.89±0.57 


0.36±0.02 
0.35±0.03 
0.34±0.02 
0.33±0.01 
0.29±0.03 


H/3 


Ob 
lb 
2b 

3g 

4g 


1.30±0.25 
1.72±0.31 
2.40±0.34 
9.51±0.47 
8.18±0.46 


I. 62±0.32 
2.14±0.38 
2.97±0.43 

II. 81±0.58 
10.15±0.57 


0.11±0.02 
0.16±0.03 
0.11±0.02 
0.24±0.01 
0.50±0.03 


[No III]A3869 


Ob 
lb 
2b 

3g 

4g 


0.92±0.22 
1.05±0.25 
1.68±0.28 
4.08±0.37 
2.00±0.36 


1.22±0.29 
1.39±0.33 
2.22±0.37 
5.41±0.49 
2.65±0.47 


0.08±0.02 
0.11±0.03 
0.08±0.01 
O.lliO.Ol 
0.13±0.02 


[O II]A3727 


Ob 
lb 
2b 
3g 

4g 


1.18±0.21 

<0.72 
2.00±0.26 
12.83±0.52 
14.30±0.49 


1.59±0.28 

<0.96 
2.69±0.35 
17.21±0.70 
19.19±0.65 


0.11±0.02 

<0.07 
O.lOibO.Ol 
0.3d±0.02 
0.94±0.04 


[Ne V]A3426 


Ob 
lb 
2b 
3g 
4g 


2.44±0.20 
2.18±0.23 
2.26±0.25 
1.24±0.32 
<0.94 


3.38±0.28 
3.03±0.32 
3.14±0.35 
1.73±0.44 
<1.30 


0.23±0.02 
0.23±0.03 
0.12±0.01 
0.04±0.01 
<0.06 


[Ne V]A3346 


Ob 
lb 
2b 
3g 
4g 


0.92±0.20 
0.79±0.23 
<0.71 
<0.96 
<0.96 


1.28±0.28 
1.10±0.33 
<1.00 
<1.35 
<1.35 


0.09±0.02 
0.08±0.03 
<0.04 
<0.03 
<0.07 



Region 0b=[7.9,13.8)" ;Region lb=[13.8, 20.7)"; Region 2b=[20.7,27.6)"; Re- 
gion 3g=[27.6, 38.63)"; Region 4g=[38. 63,49. 67)" as shown in Fig.[6l 
^ ISchleeel eFaLl flggsl: Car delli et ahl {1989) 

Dereddened fiux relative to value of [O III]A5007 in that spatial bin. 

All upper limits are at 3o" level. 



10 



dated with our data reduction (i.e. sky subtraction and 
continuum subtraction). While the former is relatively 
straightforward to estimate, the latter is more difficult. 
We estimated this systematic error through analysis of 
our science images that had been sky and profile sub- 
tracted as follows. For every pixel in these images, we 
calculated the sigma-clipped mean and standard devia- 
tion in a square region with 10 pixels on a side (after 
masking known cosmic rays and the emission lines them- 
selves). We then have an image of the local means and 
an image of the local standard deviations. For each emis- 
sion line we picked a region from these images of means 
and standard deviations that is nearby (but does not in- 
clude emission) and proceeded to take their averages. We 
then add an additional error term in quadrature with the 
Poisson errors equal to the average of the standard devi- 
ations in this region. We also include an additional error 
associated with over/undersubtracting the profiles that 
is characterized by applying the entire profile subtraction 
algorithm to regions that contain no emission lines and 
measuring fluxes in the same velocity apertures in those 
regions. 

Lastly, when comparing line fluxes we must take care to 
consider the same spatial region for each of the emission 
lines. Because our observations were taken with differ- 
ent cameras each with a different plate scale, the spatial 
region determined by an integer number of pixels for one 
of the plate scales is not generally equal to an integer 
number of pixels for others. Therefore, we insure we are 
comparing the same region by interpolating the spectrum 
onto a grid 100 times finer in spatial and velocity resolu- 
tion. We note that this is valid as our platescales, seeing, 
and dispersion render the data Nyquist sampled in all of 
the observations. We then use the spatial centroid of the 
quasar continuum light as a reference and extract along 
the same spatial intervals for all the lines. We then ap- 
ply velocity extraction windows as shown in Fig. [S) The 
velocity extraction windows are broader for the [O II] 
lines because they form a doublet and are treated as a 
single line which is therefore seen to have a larger veloc- 
ity spread. The Ha and [N II] lines have slightly smaller 
velocity windows to prevent overlap with nearby bright 
skylines. 

The flux measurements are presented in Table [31 We 
corrected fluxes for Galactic foreground reddening using 
the Schlegel, Finkbciner and Davis dust maps {Ay — 
0.292; Schlegel et a l. 1998.) and a Milky Way dust ex- 
tinction curve from iCardelli et al.l 



Milky 
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4.2. Line Ratios 

Fig. [7] presents several line ratios that gauge physical 
characteristics of the system. One such ratio, that of Ha 
to H/3 (black crosses), is a standard indicator of dust ab- 
sorption. The dotted lines show theoretical estimates of 
this ratio at a range of temperatur es spanning reasonable 
expec tations for photoionized gas (jOsterbrock fc Ferlandl 
|2006[) . The agreement between the observations and the- 
oretical prediction implies that there is little dust in the 
bridge and only modest extinction within the companion 
galaxy. 

The figure also shows how the brightness of an emis- 
sion line varies with distance from the quasar. The ap- 
parent trend is that lines coming from higher ionization 
species reach their brightest intensities (relative to other 
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Fig. 7. — Tlie spatial variation of various line ratios measured for 
the J2049-0012 system. The lines used were Ho, [O III]A5007, 
[Ne V]A3426, and [Ne III]A3869. The dotted lines denote the the- 
oretical values for the Ha/H/3 ratio for gas at Te = 5,000K (3.05) 
and at 20,000K (2.76) (Osterbrock Ferland .2006i) . The galaxy 
and bridge are consistent with having little intrinsic dust. Higher 
ionization lines are stronger relative to lower ones closer to the QSO 
(at spatial position kpc). This trend suggests that it is possible 
the quasar is responsible for the emission because the [O IIl]/H/3, 
[Ne V]/[Ne III], and [Ne V]/[0 III] ratios are all sensitive to the 
ionization state of the gas. 

lines) closer to the quasar while lower ionization species 
show the opposite trend. Thus the ratios [O III]/H/3, 
[Ne V]/[Ne III], and [Ne V]/[0 III] are decreasing with 
increasing distance from the quasar. This suggests that 
the quasar is related to the observed emission (see § 16. 2p . 

The ratio of [Ne V] to [Ne III] constrains the ioniza- 
tion parameter of the gas (which is an indicator of the 
ionization state of the gas). Indeed, the presence of 
[Ne V] emission alone requires a source of high energy 
photons or high temperature (T) gas (with an ionization 
potential exceeding 7 Rydberg, H II regions will not typi- 
cally produce this line) . The emission line ratios are also 
important in determining the source of ionization, e.g. 
[N II] /Ha and [O III]/H/3 are sensitive to the ionizing 
source's spectrum. 

4.3. Kinematics 

The kinematics of the emission constrain the driving 
gravitational forces on the gas. They also provide esti- 
mates for the masses of both the QSO host and its com- 
panion galaxy. Coupled with comparisons to simulations, 
we can also learn about the specific merger stage of the 
galaxies. To characterize the kinematics of the ionized 
gas, we calculated the flux-weighted velocity centroid and 
velocity dispersion in the [O III]A5007 line within each 
spatial row of pixels. We choose the [O III]A5007 line 
because it is the brightest observed line with the highest 
signal-to-noise ratio. These measures were calculated by 
fltting single gaussians to all velocity pixels correspond- 
ing to a given spatial row of pixels. Each fit was checked 
by eye and catastrophic errors related to non-gaussian 
line-profiles have been masked. These are due to regions 
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Fig. 8. — Kinematic characteristics for several emission lines detected for the J2049— 0012 system. The values presented on the left are 
derived from Gaussian fits to each spatial row in the 2D spectrum, {left top) The velocity centroids of each fit; km/s is defined by the 
SDSS redshift of the QSO. The arrows at the top denote the spatial extent of each region of the system. Specifically, the star denotes 
the centroid of the quasar, the circle the centroid of the galaxy and the central region is called the bridge. Note that the [Ne V] emission 
is kinematically coincident with the [O III] emission indicating that these ions likely arise in a single phase of gas. (left bottom) The 
velocity dispersions (ct) measured from the higher resolution 1200/7500 grating for the [O III]A5007 line. The dashed line is the measured 
velocity dispersion of a nearby sky line. The red dotted line represents the expected velocity dispersion as expected by the finite spatial 
resolution of the spectrum (see the text). Spatial bins that were poorly fit by a single gaussian have been masked, (right) The higher 
resolution (1200/7500 grating) spectrum of the [O III] emission. This spectrum shows finer detail than the 600/7500 grating and is used 
for comparisons with simulations and to estimate the masses of the galaxies. One can clearly see that the bridge has distinct kinematic 
properties from either galaxy, yet connects spatially with each galaxy near their velocity centroids. 



with low signal-to-noise or when there are overlapping 
velocity components (e.g., in spatial regions that include 
both the bridge and one of the galaxies) . 

The dominant source of uncertainty in these kinematic 
measures is related to systematic effects. Specifically, the 
wavelength solution based on the arc lamp images has a 
typical RJVIS residual of ~ 0.1 pixels corresponding to a 
velocity of « 10 km s~^. The systematic velocity error 
due to this effect is calculated for each line by (1) mea- 
suring the RIMS of the residuals from our wavelength so- 
lution in units of pixels (2) multiplying the RIVIS by the 
dispersion (A/pix) to get the wavelength RIVIS (3) then 
estimating the error that is introduced at each line's ob- 
served wavelength. These errors are added in quadrature 
with the fitting errors. 

We find that the kinematics of the [Ne V]A3426 line 
trace the [O III]A5007 emissior0; it is thus very likely 
that these ions are co-spatial. This is a crucial point for 
the determination of the ionization source, as it allows us 

^ It is best to compare the [Ne V]A3426 kinematics to those 
derived from the R600/7500 grating for [O III]A5007 since they 
were measured with comparable resolution. 



to demand that the same ionization mechanism produce 
each of these ions and their observed line ratio. 

The velocity centroids and dispersions (cr) of 
[O III]A5007 as measured from our highest resolution 
R1200/7500 observation are shown in Fig. [5] In our 
higher resolution spectrum, the kinematics of the gaseous 
bridge appear distinct from the rotation observed in the 
quasar's host galaxy and the companion galaxy (See 
Fig.HD. 

The lower panel reports the measured velocity disper- 
sions of the 1200/7500 measurements of the [O III] A5007 
line. The width of a nearby unresolved sky line is plot- 
ted as a dashed line to indicate the resolution limit of 
the spectrum. Since the velocity dispersions are higher 
than those of a nearby skyline we appear to have spec- 
troscopically resolved the emission of the bridge, but it 
is unclear if we have done so spatially. Specifically, we 
emphasize that if there are spatial gradients in the ve- 
locity field that are unresolved they have the effect of 
artificially increasing the velocity dispersion because the 
gradient is smeared. Thus one might misinterpret spa- 
tially unresolved velocity shear as random motions. To 
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quantify this effect we created a mock 2D spectrum witl 
the same central velocities and brightnesses as a functioi 
of position as those measured from our data. To isolati 
the effects of this spatial smearing we set the velociti 
dispersions of each spatial pixel equal to that of the un 
resolved skyline. We then smoothed this 2D spectrun 
in the spatial direction according to a spatial resolutioi 
element of 0.7" and measured the velocity dispersion ii 
the same manner as for the data. The result is repre 
sented as the red dotted curve and can be interpreted a; 
the expected effective resolution limit as a result of botl 
finite velocity resolution and astronomical seeing. Sinci 
the velocity dispersion is still higher than this curve, i 
appears we are also spatially resolving the velocity dis 
persion (as long as the gradients are not smaller than thi 
seeing) . 

The velocity of the quasar spectrum as measured fron 
AGN emission lines is generally a poor measure of th( 
syste mic velocity for its host galaxy fe.g.. lRichards et all 
[2OOI . Therefore, we have also analyzed the host galaxy'; 
extended [O III] emission. We use the measured centroid; 
on either side of the host galaxy along with the spatiai 
centroid of the quasar light to interpolate the velocity 
at the spatial center of the host galaxy. We do this by 
fitting a robust linear fit to the data. Thus we measure 
the relative radial velocities of the quasar host galaxy 
and companion galaxy to be 30 ± 30 km s~^, with the 
quasar at the higher radial velocity. 

5. PHYSICAL PROPERTIES OF THE INTERACTING PAIR 

In this section we present analysis detailing measured 
physical properties of the galaxies and the SMBH of the 
QSO involved in the merger. 

5.1. QSO Properties 

One can estimate the mass of the supermassive 
black hole Mbh fueling the quasar activity in sev- 
eral ways based on the em i ssion line charac teristics 
(jVestergaard fc Petersonll2006l : iKong et al.ll2006[ ). By fit- 
ting SED templates one may estimate the bolometric lu- 
minosity of the QSO from the continuum luminosity in 
a particular wavelength region. The quasar under study 
is in the SDSS QSO catalog of iShen et all (pOOl ) who 
used the above methods to calculate LboI = lO''^'^ erg/s 
and Mbh = lO^'^Mo where each has a typical error of 
approximately 0.4 dex. This places the quasar above the 
average for quasars at this redshift but within approxi- 
mately la of the distributions for both quantities. This 
also gives an Eddington ratio of 0.0 7 which falls within 
the normal range for quasars (iKellv e t al. 201(|). 

We utilize the cal i bratio n of iFerrarese (,20021 ) derived 
using I Bullock et al.l ()2001|) to estimate the dark matter 
halo mass of the QSO host galaxy. 
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where Vc is the circular velocity of the host galaxy at the 
point where the rotation becomes flat. 

Using the data from the R1200/7500 grisni (see Fig.©, 
we find that emission associated with the quasar host 
galaxy extends from —174.3 to -1-174.0 km s~^ (where our 
measurements at the positive velocity end are truncated 
when they become cospatial with the bridge) . These ve- 
locities however are relative to the SDSS quasar redshift 
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Fig. 9.— a BPT HBaldwin et aLlirOSl ) diagram of the flux ratios 
observed for the bridge (Ob, lb, 2b) and companion galaxy (3g, 4g) 
regions as demarcated in Fig. [6] The points are a random selection 
of galaxies from SDSS to give a sense of the normal locus of ob- 
jects. The solid curve is from Kewley et al. (2001) and separates 
regions consistent with ionization from H II regions from other 
sources of i onization . The dashed curve is a similar division found 
bv iKauffmann et al.l ll2003b ). The horizon al line segment denotes 
the se paration between LINERS and AGN iVeilleux &: Osterbrocld 
1198711 . The regions of the bridge (Ob, lb, 2b) are consistent with 
being photoionized by a hard AGN spectrum while bin 4g (located 
farthest from the quasar) is more consistent with H II regions. Bin 
3g is in a portion of the diagram that is inconsistent with ioniza- 
tion by H II regions alone and likely has an ionizing spectrum that 
includes the AGN and a stellar component. 

and not the systemic velocity of t he quasar host. Us- 
ing the methods described in M.3\ we found the offset 
to be 30 ± 30 km/s. Thus we arrive at an estimate 
of Vc > (200 ± 30)/ sin i where i is the galaxy inclina- 
tion. This is a lower limit on the actual circular veloc- 
ity because (1) we may not have aligned our slit with 
the major axis of the galaxy and (2) our spectrum does 
not appear to extend to the flat portion of the rotation 
curve. Using the la error below the measured value, we 
flnd that this gives an estimate of the host mass to be 
Mdm > 8.8 X 10" Mq. 

With IFU spectra or additional spectra at different 
slit angles, we could constrain the inclination and axis 
ratio of the host galaxy to better characterize the host 
galaxy's dynamical mass. Space-based imaging however 
would allow us to directly image the host's bulge and 
multiband imaging would allow us to better estimate its 
stellar mass. In paper II we discuss the position of the 
BH and host on Magorrian-like relations and details of 
how this and other systems like it might lead to insights 
on whether a host galaxy leads/lags its SMBH. 

5.2. Companion Galaxy Properties 
5.2.1. Star Formation Rate & Gas MetalUcity 

The observed line ratios of Ha/ [Nil] and [0III]/H/3 
define the the standard axes of the Baldwin, Phfl lips 
and Terlevich (BPT) diagram (jBaldwin et aLlflOSl see 
Fig- El) • These line ratios are used because they are sensi- 
tive to the hardness of the spectrum and the temperature 
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of the gas (which in turn is also sensitive to the hard- 
ness of the spectrum). Furthermore, these are bright 
hues in the optical passband that are closely spaced in 
wavelength such that differential reddening is generally 
negligible. In this diagram H II regions occupy a partic- 
ular locus, while AGN and LINER sources each separate 
into their own loci. The curved line denotes the division 
between emission characteristic of ionization from hot 
young s tars in H II regions from ionization from other 
sources (|Kewlev et al.ll200ir ). The dashed, curve is a sim- 
ilar division as reported by iKauffmann et al.l ()2003bl ). 
The horizontal line segment to the right of this curve 
at a value of log([0 III]A5007/7?/3) = log(3) denotes the 
division b etween the LINER and Seyfe rt line ratios as 
defined bv lVeilleux fc OsterbrocH (|1987[ ). 

The galaxy spans bins 3g and 4g. While bin 3g ap- 
pears strongly affected by an AGN-like source, bin 4g 
is closer to the locus of gas photoionized by stars. Un- 
der the hypothesis that the ionization in region 4g of 
the companion galaxy is consistent with H II regions, we 
can es timate the star formation rate using the iKennicuttI 
(|1998D calibrations for Ha: 

SFRi{MQ yr-i) = 7.9 x 10"*2L(Ha)(ergs 8-^)-^ 

= 1.61 ± 0.03 M© yr-^ (2) 

Assuming that the [Oil] emission also arises from H II 
region^, we can additionally estimate 

SFRi{MQ yr^i) = (1.4±0.4) x 10-''1L([0 II])(ergs s'^)-^ 

The two estimates are in good agreement, suggesting that 
dust is not adversely affecting the line fluxes. Of course, 
we are not sensitive to any highly obscured star formation 
and in this sense the measured star formation rate is a 
lower limit. We also note that this SFR is only measured 
over bin 4g. We also detect Ha emission for the galaxy 
within bin 3g but we expect that these line fluxes are 
influenced by the neighboring QSO's ionizing radiation. 
Therefore, we assume that the star formation in that bin 
is equal to bin 4g and apply a multiplicative correction 
factor of 2. We also need to correct for the fact that our 
slit is only 1" wide and the galaxy extends outside of the 
slit. To do this, we make use of our i?-band image and 
(1) measure the total flux of the galaxy (2) apply a 1" 
slitmask and measure the total flux which falls within 
this slit. The ratio of these two fluxes is the slit loss 
correction. We find this to be a factor of 2.5. We thus 
apply a total correction factor of 5. This of course is 
assuming that the star formation rate is constant across 
the entire galaxy and that the i?-band light follows the 
star formation. This gives an estimate of the total star 
formation rate for the galaxy of 8.05 ± 0.15Mq yr^^. 

As a caveat, we point out that the star formation in the 
companion galaxy may be dominated by its innermost re- 
gions and thus not indicative of the overall galaxy mor- 
phology. This may occur when gas is funneled by tidal 
torques of the interaction to trigger a burst of recent star 
formation. If this is the case, we may be overestimating 
the slit loss correction. For a tabulated summary of all 
SFR estimates with the varying assumptions, see Table 

Relatively weak [Oil] emission is predicted for a hard radiation 
field l|Holl2005ll . 



[31 These estimates likely bracket the range of possible 
star formation rates. 

We may also estimate the metallicity of the gas 
in bin 4g using standard line flux ratio calibra- 
tions. iPettini fc Pagell (|2004f ) outline two calibrations 
of metallicit y determination. In these units solar 
0/H is 8.66 dAllende Prieto. Lambert, fc AsDlundll2001l: 
lAsplund et al.ll2004D : 

12 + log(0/H) = 8.9 + 0.57 x iV2 = 8.64 ± 0.01 (4) 
12 + log(0/H) =8.73 - 0.32 x 037V2 = 8.49 ± 0.01(5) 

where N2 = log{[N II] /Ha} and 037V2 = 
log{([0 III]/H/3)/ ([N II]/Ha)}. These empirical rela- 
tions have scatter such that 95% of points are within 
0.38 dex of the N2 calibration and within 0.25 dex of the 
03N2 calibration. Thus the two measurements are in 
agreement within the systematics of each other and the 
solar value. Therefore, we conclude that the galaxy has 
roughly solar metallicity. 

5.2.2. Mass Estimates 

We can estimate the mass of the companion galaxy in 
the same manner as for the quasar host galaxy f i)5.ip . i.e., 
through analysis of the rotation curve. The measured 
circular velocity implies an estimate to the dark matter 
halo mass Mum > 2.1 x 10^^ Mq. This is also a lower 
limit due to uncertainties in the inclination, axis ratio, 
and position angle of the galaxy. We note, however, that 
unlike the ciuasaLjiost, the data does extend to the flat 
^8?tfefti%4fe; ¥6tation curve. 

Using the 5 SDSS colors (ugriz), one may constrain 
the spectral energy distribution of a galaxy at z < 1 to 
estimate its mass-to-light ratio. This mass-to-light ra- 
tio can then be combined with the observed luminos- 
ity to estimate its stellar mass. Mak ing use of ver- 
sion v4.1.4 of the KCORRECT code (Blanton fc Roweij 
|2007[ ). we estimate a stellar mass of 1.9 x IO^^Mq with 
a typical error of approximately 0.2 dex compared to 
other photom etric estimates of stellar m ass. Using the 
calibration of 'Conrov fc Wechsler (2009) which is very 
similar to that of Fontana ct al. (2006), we find that 
M*{z = 0.369) = 1.01 X 1O"M0. Therefore, the com- 
panion galaxy appears to be « 0.2M*. Using halo abun- 
dance matching to relate stellar mass to halo mass (e.g. 
iConrov fc Wechslej 120091 ). we estimate the dark matter 
halo mass to be approximately 6 x 10^ "'^M0. 

5.2.3. Stellar Population Models: Is the Companion a 
Post-Starburst Galaxy? 

The classic definition of a post-starburst galaxy is one 
which has strong Balmer absorption features along with 
Ca H-f K absorption lines and an absence of emission 
lines. While our galaxy does have emission lines which 
point to ongoing star formation, we believe that our 
galaxy falls into the post-starburst class for two reasons: 
(1) a significant portion of the emission may be a result 
from photoionization from the nearby quasar and (2) our 
following analysis demonstrates that the star formation 
was likely elevated in the past in order to simultaneously 
reproduce th e colors and absorp t ion lin es. 

We utilize iBruzual fc Chariot! (|2003[ ) population syn- 
thesis models to constrain the star formation history of 
this galaxy. We construct solar metallicity models using 
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Fig. 10. — Comparison of SSP models to the observed spec- 
trum. Data are shown in black, the models are shown in red, and 
masked emission line regions are shown in cyan. The observed 
strong Balmer absorption lines are characteristic of A star features 
which indicate star formation halted in a time period that is longer 
than the lifetime of O and B stars but shorter than that of A stars. 
This range is 100 Myr - 1 Gyr and we refined the age estimate with 
SSP modeling. We find that the range of reasonable models is for 
starbursts that started t=300-800 Myr ago. Older stellar popula- 
tions do not have strong enough Balmer features and are too red. 
Younger populations also underestimate the Balmer features and 
are too blue. Our best model is plotted in the first row which cor- 
responds to 800 Myr. The second and third rows show examples 
of the poor fits that result from too old and too young populations 
respectively. We also show in the bottom panel the results of using 
a star formation history determined by. Noes ke et al., £2007). This 
SFH insufficiently reproduces the spectral region between Ca K 
and H7. 



Padov a 1994 stellar evo lutionary t r acks with a ICliabri"eil 
(|2003l ) IMF. Following lYan et alJ (|2006f ). we used vari- 
ous two-component models for the burst star formation 
history. The first component is a single 7 Gyr old, in- 
stantaneous stellar population which represents an old 
passively evolving background stellar population. The 
second component is a r model burst occurring within 
the last 2 Gyr. We used a variety of r's to parameterize 
our burst of star formation [where SFR oc exp(— i/r)]: 
0.01, 0.1, 0.25, 0.5, 0.75, 0.9, 2, and 4 Gyr and spaced 
the age of the burst at 100 Myr intervals. We allowed a 
ra nge of burst amplitudes (A ; defined in the same manner 
as IKauffmann et all (|2003aD ) such that the stellar mass 
produced by the burst (evaluated at t = cx)) relative to 
the initial 7 Gyr population was between 10"'^ and 10'* 
in logarithmic steps of A log A = 0.53. 

We fit these model spectra to match our observed 
spectrum after correcting for Milky Way dust absorp- 
tion in the same manner we treated our emission flux 
measurements (see § 14. ip . Due to the LRIS dichroic 
gap we restricted ourselves to two rest wavelength ranges 



3700-4000 A and 4250-4900 A. We also masked ah pix- 
els within 400 km/s of an emission line since our models 
do not include nebular emission. We allowed the model 
spectrum to scale up or down by a single constant fac- 
tor and allowed for an intrinsic dust corre ction. The dust 
correction was implemented with a . Chariot fc Falll (|2000f ) 
law parameterized by ry: 
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where -Fx ohs is the observed fiux and -Fx i„f is the intrinsic 
flux. 

We applied the additional criterion that the SDSS w— g, 

— r, r — i, and i — z colors of the models (including ap- 
plying the Milky Way dust absorption and fitted intrinsic 
dust) are within 0.5 mag of the observed colors. 

Altogether, we created 2016 models. We find that no 
individual model provides a good (x^ — 1) fit to all of 
the spectral features. However the models that simulta- 
neously best reproduce the observed colors and the line 
features (both the Balmer and Ca H-l-K) have ages for 
the burst component that fall in the range of 300 - 800 
Myr, star formation histories parameterized by 0.25 Gyr 
< T < 0.9 Gyr, and have a burst mass fraction (rela- 
tive to the old passively evolving component) that fall in 
the range 0.009-0.1 (see Fig. ITO)) . However we note that 
among those models there are a handful which have a 
mass fraction as high as 4. All models have dust red- 
dening characterized by Ty < 1.01 with very few hav- 
ing Ty > 0.7. The best model has Ty — 0, consistent 
with the observed Ha/H/3 ratio (Figure[7]). Older models 
are too red, while younger models are too blue. Higher 
burst fractions underestimate Ca H+K and g-band ab- 
sorption features while lower burst fractions have insuffi- 
cient Balmer features or too red colors. In summary, we 
find the data are best described by a substantial burst 
that occurred 800 Myr ago and that models with t < 300 
Myr and t > 800 Myr are disfavored. 

A reasonable question to ask is whether the companion 
galaxy is consistent with a normal, non-burst galaxy. To 
do this we make use o f the star formation histories de- 
rived bv lNoeske et al.l ()2007[) for star- forming galaxies at 
the redshift of our system. They present parameters of r 
model star formation histories as a function of baryonic 
mass. We have an estimate of the stellar mass of the 
galaxy (see above) which we assume to be the baryonic 
mass of the galaxy; corrections to include gas up to rea- 
sonable gas fractions (< 50%) have little to no effect on 
the following results. Thus, we find the star formation 
parameters to be r = 25.06 Gyr and after converting 
from redshift intervals to ages, we find the beginning of 
the model be 5.13 Gyr prior to our current observation 
of the galaxy. Thus this model is very nearly a constant 
star formation rate model. We fit this model to the data 
allowing for a normalization factor and dust reddening 
as we did for the burst models. We present the results in 
the bottom panel of Fig. [TUl One can see that this model 
insufficiently produces the spectral features between Ca 
K and H7. Additionally, the model's large amount of 
ongoing star formation produces too blue of an SED and 
thus requires a rather large column of dust to redden the 
spectrum (ry = 1.27). Assuming the same amount of 
reddening corresponds to the H II regions that we use 
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to estimate the SFR, we find two results. Firstly, this 
measure of reddening implies a higher Ha/H/3 ratio than 
observed, but only by > la. Secondly, the reddening 
correction would increase our estimate of the SFR by 
a factor of 3.7 giving a value greater than 20 Afo/yr. 
Comparison with the ob servational results o f the specific 
star formation sequence (jNoeske et al.ll2007D . we find this 
star formation rate lies well above the 84th percentile for 
galaxies at this stellar mass. These pieces of evidence 
argue against this scenario, but we cannot rule it out 
altogether. 

The relatively high level of ongoing star formation may 
also be an indicator of a previous phase of elevated star 
formation. S FHs for galaxy me rgers are not (5-functions. 
For example, iCox et all (|2006f) show that while the SF 
activity reaches a peak after first passage it can remain 
elevated until the merger ultimately coalesces. 

6. ANALYSIS OF THE BRIDGE 
6.1. Density & Mass 

We can make an order of magnitude estimate of the 
column density along the line of sight through the bridge 
by using the observed luminosity of the bridge in recom- 
bination lines combined with a simple volume model for 
the bridge (See Fig. [5]). We note that the column den- 
sity we are about to estimate is not the column density 
of material that the quasar radiation encounters along 
the bridge, but is the projected column on the sky from 
our perspective. We denote this column density as Nh.i_ 
because if the bridge is in the plane of the sky then this 
column density is perpendicular to the column density 
that the quasar's radiation encounters. 

We start with the expression for the luminosity of the 
Ha line produced by recombinations: 

Lua neUpolllv fhvna (7) 

where Lho is the luminosity of Ha in a given spatial 
region, Ue is the number density of electrons, Up is the 

number density of protons, a^^^ is the Ha effective re- 
combination coefficient, V is the volume the gas fills, / is 
the filling factor of the gas in that volume (observation- 
ally unconstrained but by definition < 1), /i is Planck's 
constant, and VHa is the frequency of the Ha transition. 
We assume that the bridge is a cylinder with diameter 
T) that extends a radial distance ri from the QSO to a 
larger radial distance r2 (we note that r'i,r2 3> 2?; see 
Fig. El). In this case, V {r2 - ri). 

Assuming that the hydrogen is fully ionized, we take 
ne = np = uh, where nn is the number density of hy- 
drogen in all ionization states. This is a reasonable as- 
sumption given the observed line ratiofQ- 

Another piece of evidence that the gas is mostly ionized is 
from the observation that the emission extends along a long spa- 
tial extent. If the quasar is responsible for the ionization ( i|6.2ll the 
gas must therefore be optically thin to ionizing radiation along 
a pathlength of order r2. The optical depth to ionizing pho- 
tons for neutral hydrogen reaches unity for column densities of 
~ 10^^ cm~'^. Since our column density is measured along a path- 
length of I? <5C r2 , as long as we measure a total hydrogen column 
density (H I and H II) of > 10^^ cm~^ then this is a safe assump- 
tion. If there is more than enough gas to be optically thick if not 
ionized and the gas is not optically thick, then the gas must be 
ionized. 



Manipulating eqn. [7] and substituting for our expres- 
sion of V we define: 

_ -fl/2 / Lho -,„21 -2 

uj = J ' uhV — , / ^ w 10 cm . 

V ('"2 - ri) a^^hvHa 

This quantity consists only of observed quantities and 
atomic constants. The order of magnitude estimate is 
made using the observed flux of Ha for bin lb which we 
round to 4 x 10~^^ erg/s/cm^. Using the observations 
from bin lb, we put a limit on the volume density of the 
gas by noting that: 

Taking the limit to V placed by our Kast observations 
(< 10 kpc), we derive nn 0.1 /~^/^ cm~^. This is 
consistent with volume densities estimated for similar 
gaseous bridge structures in galaxy merger simulations 
by Weniger et al. (2009), where they implemented a mul- 
tiphase ISM code in a galaxy merger simulation. 
We may now estimate the column density, 

^ # of Hydrogen atoms ^ fnnV 

^'^ Surface Area those atoms fill V (r2 — ri) 
^nnfV « « /1/2 1021 cm-2, (10) 

This value is consistent with that of 21 cm neutral hy- 
dr ogen observations of t idal tails of the Antennae galaxy 
bv lHibbard eTall (l200l) . If that gas were fully ionized, 
then it would give rise to a similar column density of H+ . 
Lastly, we estimate the hydrogen mass 

Mh - nnVfrnp ^ V {r2 ~ n) ujf^'^rup 

This is relatively modest, but we note that T) could be 
as large as 10 kpc although / could also be significantly 
smaller than unity. 

6.2. Source of Ionization 

We observe the bridge in emission lines from recombi- 
nation of Hydrogen and forbidden lines of high ioniza- 
tion species. This, of course, requires a physical source 
of ionization and heating of the gas. In this section, we 
compare our observed line ratios against models to deter- 
mine the source of ionization. After ruling out collisional 
ionization and photoionization from fast radiative shocks 
and H II regions, we conclude that the quasar is shining 
on the bridge and photoionizing the gas and test this 
hypothesis using cloudy models. 

6.2.1. Shock Ionization 

Shocks are common occurrences in galaxy mergers and 
also result from the interaction of quasar jets with the 
interstellar and cir cumgalactic media of the galaxy (e.g. 
iRosario et al.l2010t ) . If those shocks are fast enough, they 
could ionize gas along the bridge. There are two ways in 
which shocks can lead to ionization: (1) collisional ion- 
ization of the shocked gas or (2) radiative shocks that 
photoionize the pre-shock gas with the radiation emitted 
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Fig. 11. — Comparis on of the predict ed line ratios for fast, ra- 
diative shocks (curves; lAIIen et al.ll200H) with the observed values 
in the bridge (bins Ob, lb, and 2b). These predictions are for fast 
radiative shocks+precursor in gas with a density of 1 cm~^ and so- 
lar metallicity. The range spanned by the predictions is relatively 
insensitive to these choices. The bridge's emission is inconsistent 
with shocks even if we allow those shocks to be unreasonably fast 
(up to 1000 km/s). 

by the post-shock gas. We begin by addressing the first 
possibility: whether gas in the bridge has been coUision- 
ally ionized. 

If the gas were cohisionally ionized, then emission from 
high ionization potential states such as Ne V, would re- 
quire a very large electron temperature. We can estimate 
Tg using the [O III] transitions as these form a temper- 
ature diagnost ic in the dens ity regime Ue < 10^ cm~'^ 
(fOsterb rock fc Ferlandll2006l ): 



JA4959 + JA5007 7.9 exp(3.29 X lO^Te) 



JA4363 



1 + 4.5 X lO-'^njTt 



1/2- 



(12) 



This diagnostic is valid as long as the density is below 
the critical density of the [O III] transitions. If the den- 
sity were higher than this critical density, the flux would 
be attenuated dramatically since at these densities col- 
lisional deexcitation dominates over radiative deexcita- 
tion. Such high density, however, would imply much 
lower [O III]/H/3 ratios than observed. Such a density is 
also many orders of magnitude higher than normal ISM 
densities and is highly unlikely. The conservative upper 
limit to the [O III]A4363 line flux sets an upper limit to 
the gas temperature of Te < 7 x 10'' K. This is over 4 
times lower than the characteristic temperature of Ne V 
(3 X 10^ K; Szentgyorgyi ct al. 2000). The detection of 
[Ne V]AA3426,3346 combined with the temperature up- 
per limit from [O III] rules out the hypothesis that the 
primary ionization source is collisional ionization. 

Radiative shocks, via photoionization, can produce 
high ionization lines in conditions with lower electron 
temperatures. In such shocks, the post-shock gas reaches 
high temperatures and its emission photoionizes the pre- 
shock gas. This produces a gas that is ionized by a 
relatively hard spectrum. However, the observed emis- 
sion line ratios in the bridge are inco nsistent with those 
given by the diagnostic diagrams of lAllen et ahl ([20081 ) 
for fast radiative shock -|-precursor models that span 



shock velocities of 100-1000 km/s, magnetic paramters 

B/Vn = 10"'' - 10 ^iG cm^/^ n = 0.01 - 1000 cm^^, 
and a variety of abundance ratios. These models take 
into account both emission of the shocked material and 
emission from the pre-shock gas. The line ratios that 
are clearly discrepant with our observations include the 
[O II]A3727/[0 III]A5007 and [Ne V] A3426/[Ne III] 
A3869 ratios. Most striking is that all of these models 
predict [Ne V]A3426 fainter than [Ne III] A3869, in direct 
contradiction with our observations of the bridge (see 
Fig.n]). 

Lastly, we stress that the bridge emission extends along 
the entire length of the > 38h^2 kpc bridge connecting 
the two galaxies. It could be difficult to simultaneously 
generate shocks across such a large spatial extent. Fur- 
thermore, if we assume that the beginning of this inter- 
action (and hence when the shocking began) is consistent 
with the age that we estimate from the stellar popula- 
tion modelling, this could require shocks to be sustained 
over a very long spatial extent and for a very long time 
(300-800 Myr). 

Using all of the above reasoning, we rule out shocks as 
a viable mechanism for the ionization of the bridge. 

6.2.2. H II Regions 

Studies of merging systems often reveal evidence for 
new star formation, and therefore elevated populations 
of young and massive stars. These star-forming regions 
are expected to be surrounded by H II regions of ionized 
gas which could produce the line emission observed in the 
bridge. We have tested this hypothesis in the following 
manner. 

Looking at the BPT diagram (Fig.lH]) as first presented 
in § 15.2.11 we see that the spatial bins that contain gas 
in the bridge (Ob, lb, 2b), show the bridge emission line 
ratios are consistent with gas ionized by an AGN spec- 
trum. The gas in these bins also emits [Ne V] lines; com- 
mon tracers of AGN photoionization (jAbel fc Satyapall 
l2008f l. In fact the detection of Ne V strongly rules out 
H II regions as the dominant ionization source because 
stars lack a hard enough spectrum to effectively ionize 
Ne to this high of an ionization state. 

Looking at the two bins which cover the galaxy (3g and 
4g), the emission characteristics vary dramatically across 
the two regions. While bin 4g (the far-side of the galaxy) 
has a strong contribution from H II regions, bin 3g (the 
near side of the galaxy) falls somewhat in between the 
regions of AGN and H II region ionization. This abrupt 
change in emission characteristics is consistent with a 
model where the quasar is photoionizing the gas in the 
bridge and the near side of the galaxy but where the 
galaxy effectively shields its own far side from the radia- 
tion. The abrupt change of the emission characteristics 
of bin 3g from those of the bridge is likely the result of the 
gas changing from the optically thin conditions of bins 
Ob, lb, and 2b to the optically thick conditions of the 
disk of the galaxy in bin 3g. Bin 3g, also likely has some 
contribution from H II regions, complicating its analysis. 

6.2.3. Quasar Ionization Modelling & Bridge Metallicity 

In the previous sections we ruled out shocks and stel- 
lar photoionization as the source of ionization for the 
bridge material. We also inferred, from investigation of 
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Fig. 12. — Plots revealing the allowed parameter space in ioniza- 
tion parameter U and metallicity [M/H] for the emission observed 
in the bridge (bins lb and 2b). The grey region in each plot is ex- 
cluded because these models predict significant emission from lines 
that are not detected at that level (e.g., Fe VII, O VI, O II). The 
colored areas in each figure represent the la contours for a series 
of line ratios as labelled. The black and yellow contours correspond 
respectively to the la and 2a confidence intervals for the combined 
contours of all 4 line ratios. Note that the entire region which is 
not grey is allowed at la confidence for the [Ne V]A3462/H/3 ratio 
for bins Ob and lb. Table [4] lists the best fit values. 

the BPT diagram, that the hne ratios were hkely the re- 
suh of photoionization by a quasar. In this section, we 
test this hypothesis by comparing our observations with 
photoionization models. 

Emission hne ratios can inform one about the physical 
conditions of a photoionized gas as well as the properties 
of the incident ionizing radiation. We utili ze CLOUDY 
version 08.00 photoionization simulations (F erland et alj 
[1998) to deter mine if photoionization by a quasar's ra- 
diation field can reproduce our observed line ratios and 
then we infer physical conditions from these ratios. 

Our observations include a variety of line ratios 
for this purpose. We focus on lines with relatively 
high signal to noise and ratios that depend sensitively 
on ionization state and metallicity: [O III]A5007/H/3, 
O III]A5007/[Ne V]A3426, [Ne V]A3426/H/3, and 
Ne V]A3426/[Ne III]A3869. By studying hne ratios in- 
stead of absolute luminosities, the results are insenstiive 
to the detailed geometry of the system. 

Our models are plane-parallel slabs of gas illuminated 



by an ionizing spectrum. We parameterize the shape of 
the ionizing continuum as a power law with a spectral in- 
dex of a = —0.4 from 0.01 to 20 Rydberg (consistent with 
results from modelling the QSO SED as described in Pa- 
per II). The results are not very sensitive to changes in the 
spectral range the power law covers. We included the ap- 
propriate CMB as well as cosmic rays, but these sources 
of ionization have little consequence for our calculations. 
We ran models with metallicities of —2.5 <[Fc/H]< 1.0 . 
We note that our implementation of metallicity is just a 
scaling of all elements equally relative to solar and hence 
[Fe/H] can be interpreted as a generic [M/H], [O/H], or 
any [X/H]. We varied the strength of the ionizing ra- 
diation by varying the ionization parameter U which is 
defined as the ratio of the ionizing photon density to hy- 
drogen density: 



log U = log 



S 



47rr^ 



cn 



(13) 



where r is the distance from the ionizing source, n is the 
number density of the gas, and S is the ionizing photon 
luminosity of the source. We set the hydrogen density to 
be lOcm""^ (however the particular value was not rele- 
vant because we varied the ionization parameter U, and 
optically thin gas is homologous in U). We chose the gas 
stopping column density to be such that the gas was op- 
tically thin {Nhi = 16.8). We have adopted the optically 
thin regime because we see ionization extending through- 
out the entire 38 kpc extent of the bridge. Optically 
thick models would need to be clumpy and only partially 
covering the source. We discuss futher the implications 
of clumpy models in Paper II. 

We compared our line ratios to the models using a 
statistic which requires careful estimation of the errors on 
those line ratios. We refined our error estimates by in- 
cluding an additional 10 percent error on lines if they are 
part of a line ratio that spanned different cameras to ac- 
count for relative fluxing issues. The models do not take 
into account any dust that may be a result of an inter- 
vening screen of material between the region of emission 
and the observer. To characterize the error associated 
with this, we turned to our measured line ratios. Specifi- 
cally, we characterized our errors from reddening through 
1000 Monte Carlo realizations of the observed Ha/II/3 ra- 
tio. We adopted the observed ratio with a normally dis- 
tributed error to determine the intrinsic reddening and 
then corrected the line ratios. This created 1000 "actual" 
line ratios whose standard deviation we used to estimate 
the error due to intrinsic reddening. The log relative 
errors on ([O III]A5007/H;3, [Ne V]A3426/[Ne III]A3869, 
[Ne V]A3426/H;3, [O III]A5007/[Ne V]A3426) were (-1.4,- 
0.9,-0.12,-0.62) for Ob, (-1.8,-1.5,-0.8,-0.7) for bin lb, and 
(-4.1,-3.7,-3.0,-3.05) for bin 2b. The relative errors are 
much lower for bin 2b, because the observation places it 
slightly below the theoretical value. 

We first ruled out any models that predicted detectable 
fluxes (30% of the luminosity of [O III]A5007) for high 
ionization lines (e.g. Fe VII 6087, O VI 5291, etc.) or low 
ionization lines (e.g. O II) which are not detected at that 
level in the bridge. This rules out low and high ionization 
paramters (logf/ < -3;logC/ > -2) where [OIII] 5007 
is predicted to be weak relative to these lines. This re- 
sulted in the excluded region denoted grey in Fig. [121 



18 



TABLE 4 

Measured & Derived Properties 



Property 



Estimated Value 



Merger 



Section 



2QSO 

observed Av 
Projected Distance 
Merger Stage 



0.369 ± 0.001 
30 ± 30 km/s 
38 

Between First and Second Passage 



M. 
3J 

m 



Mbh 






Shen et al. (2008) 


LboI 


1045.7 gj.g/s 


5.1 


Shen et al. (2008) 


Eddington Ratio 


0.07 


5.1 


Shen et al. (2008') 


Host DM Mass 


> 8.8 X 1O"M0 






Companion Galaxy 



Bridge 



Concentration 


>2.61 





SFR4 


1.61 ± 0.03Mq yr-l 


^.'l.i 


Assuming SFR3=SFR4; SFR 


3.22 ± 0.06M© yr-l 


5.2.1 


Assuming R-band Light traces SF; SFR 


8.05 ± 0.15M© yr-l 


5.2.1 


12 + log(0/H) 


8.64 ±0.01 ±0.2 


5.2.1 


Mdm 


> 2.1 X 10"Mq 


5.2.2 




1.9 X IQIOjUq 


^.'l.'l 


SHAM Mdm 


6 X 10" M0 


5.2.2 


SSP Age 


300-800 Myr 


5.2.3 


Burst Mass Fraction 


1-100% 


5.2.3 



Kinematics 


See Fig. [8] 

< W.2h^^ kpc 


[431 


V 


[331 




leu 


Nh,± 


« 1021 cm-2 




Mh 






Source of Ionization 


QSO Photoionization 


(Ml 


Bin Ob log U 


r, 40+0.03 


|6.2.3| 


Bin Ob [M/H] 


'J- J -0.16 


|6.2.3| 


Bin Ob xl 


0.59 


|6.2.3| 


Bin lb log U 


-2.50l°:g3 


|6.2.3| 


Bin lb [M/H] 


-0.36l°/J 


|6.2.3| 


Bin lb xl 


0.47 


|6.2.3| 


Bin 2b log U 




|6.2.3| 


Bin 2b [M/H] 


-0 23+° " 
'-'■^•^-0.16 


|6.2.3| 


Bin 2b xl 


1.32 


|6.2.3| 



With well-characterized errors and the suite of models 
in hand, we next constructed contour surfaces as a 
function of both U and [M/H]. Here we see that the 
[Ne V]A3426/[Ne III]A3869 constrains \ogU to a narrow 
region. We see that our best constraint on [M/H] comes 
from [O III]A5007/H/3. The intersection of these two al- 
low a that region is double valued. The overall and 
the [O III]A5007/[Ne V] ratio, however, break this degen- 
eracy. Due to lower signal-to- noise, the [Ne V]A3426/H/3 
ratio offers little additional constraint. The combined 
constraint lies in the same region for all three bins and 
allows us to estimate both [M/H] and \ogU with rela- 
tively high precision (summarized in Table . 

The best fit models for bin Ob, bin lb and bin 2b 
had Xr of 0.59, 0.47 and 1.32 respectively and predict 
[O III]A4363/ [O III]A5007=0.03 which is weU below our 
detection limit. Therefore these models reproduce the 
observed [Ne V] emission at a low electron temperature 



(Te ^ 1.8 X 10** K). There are additional systematic er- 
rors associated with the simple geometry in cloudy, the 
presumption of equilibrium, and uncertainty in the ion- 
ization spectral shape (e.g. variations in the spectral in- 
dex of the power law a). We characterize a systematic 
error based on variations of the power law index of a few 
tenths to produce an additional uncertainty on log U of 
w 0.05 dex and on [M/H] of w 0.06 dex. 

Lastly, we need to determine if the derived values for 
the material make physical sense. Given th a.t a quasar 
can o utput 10^^ ionizing photons per second ()TadhunteH 
Il996f ) and the bridge is approximately 20 kpc away, this 
yields a possible U parameter of \ogU — —1.16 — logn. 
Noting that the QSO could always be partially obscured 
or non-isotropic, the quasar is a feasible source for any 
density greater than logn « —1.5 cm~^. We also point 
out that our estimates for U do not appear to fall off as 
1/r^ from the quasar. In Paper II, we discuss various 
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possible interpretations of this observation. 

7. DISCUSSION 

In this section, we first synthesize the primary re- 
sults of the preceeding sections (summarized in Tabled) 
and then consider the implications for research in galaxy 
mergers and evolution. 

The J2049— 0012 system is a pair of interacting galax- 
ies. Its hallmark is an emission-line bridge that con- 
nects the two galaxies, observed in recombination and 
forbidden line-emission. The observation of the bridge 
requires that the galaxies (now with a projected separa- 
tion of 38 kpc) have had at least one tidal interac- 
tion. One of the galaxies is currently in a quasar phase 
{LboI = 10^^-^ erg/s) and its ionizing radiation has pho- 
toionized the gas in the bridge, producing the observed 
line-emission. Using a simplistic volume model for the 
gas in the bridge, we constrained its surface and volume 
densities and provided an estimate of the total gas mass 
(-IQSMo). 

Despite its neighbor undergoing a violent quasar phase, 
the companion galaxy shows no detectable AGN activ- 
ity. Its photometry and spectrum place it within or 
near the galaxy 'green valley' and also suggest a de- 
clining starburst stellar population (with burst age of 
300 — 800 Myr). The companion galaxy has a cur- 
rent SFR of ^ 6 Mq yr^^ and a gas-phase metallic- 
ity consistent with solar. This metallicity is higher by 
« 0.4 dex than our estimate for the bridge, as derived 
from Cloudy pliotoionization modelling. 

We made several mass estimates and constraints for 
both the quasar host and companion galaxy. In the con- 
text of galaxy mergers, the mass-ratio of the galaxies 
is a fundamental quantity because it strongly influences 
the dynamics of the interaction. Our system has several 
pieces of evidence which suggest it has a mass ratio of « 
1:4 and perhaps somewhat smaller. Firstly, examination 
of the dynamical measurements shows that the compan- 
ion galaxy's rotation curve flattens at a velocity of 110 
km s~^ while the QSO host's extends to a velocity of 200 
km s~^ and shows no signs of flattening. Unless the com- 
panion galaxy is seen relatively close to face on, it likely 
has a significantly lower DM mass than the QSO host. 
Additionally, our estimates for the stellar masses of the 
two systems is roughly 1:4. For the companion galaxy, 
we analyzed the SDSS photometry to estimate a stellar 
mass of 1.9 x IQ^'^Mq. For the galaxy hosting the quasar, 
we cannot directly measure the stellar light but instead 
estimate its stellar mass using standard SMBH-host re- 
lations. Taking t he estim a.ted SMBH mass of lO^ -'^Mg^ 
(|Shen et al.ll2008D and the iMagorrian et all (|1998[ ) rela- 
tion Mbh ~ 0.006 Mhost,*, we estimate 8 x 10^^ Mq 
for the stellar mass of the bulge component. This gives 
a stellar mass ratio of «1:4. If the quasar host galaxy 
has a stellar mass in a disk component then this ratio 
shrinks even smaller. This result also hinges on whether 
the black hole currently lies on the Magorrian relation. 
We conclude, therefore, that the system has a mass ratio 
of approximately 1:4 with an uncertainty on the order of 
a factor of 2. Therefore, the event most likely either just 
satisfies the major-merger (1:4) definition or lies some- 
what below the canonical threshold. 

Our observational analysis of the J2049— 0012 system 



provides a comprehensive and quantitative description 
of an interacting pair of galaxies at 2 ~ 0.4. In turn, 
it affords the opportunity to directly test the standard 
paradigm for galax y mergers as modeled by numeri- 
cal si mulations (e.g. 'Chilin garian et 311120101 : ICox et aD 
|2006[) . The current theoretical understanding of galaxy 
interactions is that galaxy mergers with mass ratios 
larger than « 1:10 generally follow a basic sequence of 
events with two possible phases of elevated star formation 
and AGN activity: one right after first passage and the 
other during final coalescence. The absolute and relative 
strengths of these two phases depends on the galaxy mor- 
phology (e.g. whether the galaxies have gas, the presence 
of a bulge) and details of the star-formation prescription. 
In such a paradigm, the combination of the detection of 
a tidal bridge and the current projected separation of 
38h^2 kpc place a strong constraint on the merger stage. 
Numerical simulations of major mer gers with a variety of 
orbits (e.g. iChilingarian et al]|2010D show that after first 
passage the galaxies move to a distance of approximately 
100 kpc where they reach their point of maximal separa- 
tion (first apogee). Subsequent passages and separations 
do not achieve separations as large as 38 kpc. This result 
is nearly independent of gas fraction ()Lotz et al.ll2010al) 
and mass ratio as lon g as the mass ratio is larger than 
1:9 ()Lotz et al.ll2010bD . We conclude, therefore, that the 
J2049— 0012 merger is after first passage (the bridge re- 
quires at least one tidal interaction) and before second 
passage. 

As described above, the current paradigm predicts that 
gas-rich galaxies will undergo a burst of star-formation 
that initiates during or just after first passage (e.g. 
iMihos fc Hernqulstl 119961 ) . This is also the only time 
aside from coalescence where gas inflow is predicted to 
be sufficiently strong to inspire strong AG N activity 
(jSpringel et al.l l2005at iHopkins et all l20Q8aD . Regard- 
ing the J2049— 0012 merger, we identify examples of 
each: one galaxy hosts a bright quasar and the com- 
panion galaxy shows signatures of a starburst stellar 
population (Figure ITOl) . Consistent with theory of in- 
teractions causing tidal stripping and strong gaseous in- 
flows which lead to large central star bursts, the com- 
panion galaxy is observed to be consistent with being 
centrally concentrated (we measure a concentration of 
at least 2.61, limited by the image's spatial resolution). 
Such compact companions to quasars have been observed 
around other more n earby QSOs (Canalizo &: Stocktonl 
120011: lStock"tonlll982| ). With standard stellar population 
modeling, we found the burst occurred 300 — 800 Myr 
ago. This timescale is significantly shorter than the pre- 
dicted period between first and second passage for galax- 
ies with our allowe d range of mass ratios for standard or- 
bits (Ki 2 Gyr; e.g. iLotz et al.ll2010bl : IChilingarian et al.l 
I20T0I) . 

In merger simulations, the level of enhanced star for- 
mation during the first passage starb urst is sensi tive to 
the star formation prescription (see iCox et ahl [2006). 
Therefore characterizing this phase may allow greater in- 
sights on sub-grid models. To this end, we would like to 
constrain the observed mass fraction of the burst. Unfor- 
tunately, various degeneracies limit us to only constrain 
the burst to consist of between 1 — 100% of the stellar 
mass prior to the merger. 
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TABLE 5 
Comparison with Theory 



Theoretical Prediction 



Reference 



Observation 



Consistent? 



Galaxy internal parameters 
determine their evolution in a merger 

First passage inspires a starburst 



Mergers may lead to quasars 

This quasar phase typically 
occurs during final coalescence 

Starburst & tidal effects may 
lead to a concentrated morphology 



2 Gyr timescale between first 
and second passage 

Kinematics near apogee have 
small relative velocity 



Mihos fc Hernauist (1996) 

Mihos fc Hernguist (1996) 

Hopkins et al. (2008a) 
Hopkins et al. (2008a) 



Stockton (1982) 



Chilingarian et al. (20101 
i.e. GalMer 



One galaxy is a quasar, 
one has no AGN activity 

Starburst spectrum observed 
after first passage, before second 

One galaxy observed as a quasar 

Merger between first 
and second passage 

Imaging is mildly concentrated, 
but inconclusive due to 
insufficient spatial resolution 

t = 300 - 800 Myr after first 
passage, before second passage 

Av = 30± 30 km/s 



Y 



Y 
7 



Another test of the merger paradigm combines the tim- 
ing of the starburst and the observed configuration of the 
merger. The current projected separation of 38 kpc 
combined with the timescale of 500 Myr requires that the 
galaxies have been moving apart from each other at an 
average speed of approximately 75 km s^^. This is con- 
sistent with the observed timescales and distances seen 
in GalMer merger simulations (Chilingarian et al. 2010). 
The observed offset in radial velocities of the two galaxies 
also follows the standard picture. The centers of mass for 
the two galaxies (as traced by the bright emission lines) 
are moving at a relative radial velocity of 30 ± 30km s~^. 
This is consistent with the system being close to apogee, 
as predicted for a system with an ag e of i ^ 500 Myr post 
first passage (e.g. iLotz et al.l[20T0bl ). 

Lastly, merger simulations also offer robust predictions 
on the physical characteristics of the tidal bridge. In the 
following few paragraphs, we discuss how measurements 
of the bridge properties test the merger paradigm and 
also present a means for gaining greater understanding 
into specific mergers. 

Firstly, we note that the metallicity of the bridge, 
which we find to be a few tenths dex lower 
(A[M/H]=— 0.36 ± 0.23) than the companion galaxy, is 
consistent with having been tidally stripped from the 
(presumed) lower metallicity outskirts of that galaxy. 
This is the standard result of the merger paradig m first 
studied in detail by iToomre &: Toomrd (|l972bl ). The 
outer parts of these extended tidal features arise from 
the outskirts of galactic disks because the tidal forces 
are greatest in these regions. 

Due to the unique emission line bridge of our sys- 
tem, we are able to measure the gas surface density 
at a point midway between the two interacting galaxies 
(i j6.ip . Wishing to use this measurement to gain more 
insight on the interacting pair of galaxies, we turn to 
numerical simulations. As this is not a standard observ- 
able and lacks sufficient literature detailing its evolution, 
we have explored the GalMer databas e of galaxy mergers 
to characterize the bridge's evolution (jChilingarian et all 



l2010r ). Using a suite of 1:1 and 1:10 galaxy merger simu- 
lations with a variety of spiral morphologies (1:1 Sa-Sa, 
Sa-Sb, Sa-Sd, Sb-Sb, Sb-Sd, Sd-Sd ; 1:10 SO-Sa, SO-Sb, 
SO-Sd ), we measured the gas surface density in the cor- 
responding regions of these galaxy mergers at 64 view- 
ing angles evenly spaced in solid angle at all applicable 
timesteps (where the galaxies are separated by at least 
20 projected kpc and have undergone at least one tidal 
interaction). The 1:1 mergers have 12 different orbits 
with each orbit also having 4 different relative orienta- 
tions of the disks. The 1:10 mergers are limited by the 
availability of the models to only have 6 orbits with only 
1 generic disk orientation each. 

We find that for 1:1 mergers, the bridge forms after first 
passage and that its surface density reduces roughly ex- 
ponentially with time (approximately oc exp[— i/r] with 
T «100 Myr and the normalization correlating with the 
gas mass of the two galaxies) as the material rains back 
onto the two separating galaxies. In the vast majority of 
cases, the gas density falls far below the measured surface 
density of our bridgeQ before reaching second closest ap- 
proach. These findings indicate standard merger scenar- 
ios of galaxies with high gas fraction can reproduce our 
observations provided we have caught the system before 
maximal separation. 

For 1:10 mergers we find that the gaseous bridge does 
not reach sufficiently high surface densities in the bridge 
for any of the simulations we examined. We do note 
that in our 1:10 simulations the larger galaxy does not 
contain any gas (we are limited by availability of sim- 
ulations), which is likely a poor match to our system. 
This is unlikely to affect our overall results because in 
such a high mass ratio merger the contribution of the 
more massive galaxy to the extended tidal features will 
be minimal because the tidal force it feels is relatively 
small. In such a case, the smaller galaxy is much more 

We note that while our estimate of the column density is 
dependent on the filling factor (oc f^^'^), the dissipation is rapid 
enough that our result is relatively robust. We use our value of the 
upper limit corresponding to / = 1. 
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greatly tidally perturbed, resulting in extended tidal fea- 
tures but at much smaller surface density than observed. 
If the mass ratio of the J2049— 0012 system were as small 
as 1:10, these models would be ruled out. 

The kinematics of the bridge further enlighten us about 
the system (See Fig. [8|). The velocity centroids of the 
bridge as a function of spatial position show that the side 
of the bridge closer to the companion galaxy has a higher 
radial velocity than the side closer to the quasar. If the 
gas is freely flowing towards the centers of mass of the 
two galaxies, this suggests that the companion galaxy is 
background to the quasar. This argument is confirmed 
by an examination of the GalMer merger simulations. 
The measured velocity dispersion of 50 km/s indicates 
an intrinsic velocity dispersion of 40 km/s after we re- 
move instrumental effects. This is much higher than that 
found in 21 cm maps of tidal debris around local inter- 
acting galaxie s which have typical velocity d i spers ions of 
10 -20 km/s (|Hibbard fc van Gorkomlli993L Il996h . One 
possible explanation is suggested by the clumpy nature 
of the emission: portions of the bridge may have become 
self-gravitating. Such structures are observed in the tidal 
features of local ga laxies and can have large velo city gra- 
dients > 40 km/s (|Weilbacher et al.ll2003l l2002l ). These 
structures are also a pred iction of merger simulations 
(|Barnes &: Hernguistl 119921 ) . It is possible that one or 
several of these structures could add to the velocity dis- 
persion. One immediate challenge to this possibility is 
the lack of obvious star formation associated with the 
bridge. However, if the bridge is fully ionized, it is possi- 
ble that the large ionization fraction of the gas has sup- 
pressed star formation. This is supported by the recent 
observational evidence that there is a neutral gas col- 
umn density threshold (IorNh t > 20.6 cm~^) for star 
formation in tidal material (Ma vbhate et al. 2007). An- 
other possible explanation is that dispersion has been 
elevated by some sort of turbulence possibly associated 
with shocks that are thought to occur throughout the 
course of the merger as a result of the bulk motions of 
the galaxies. 

Given that we observe a bridge of material that we 
believe to be photoionized by the quasar, one might ask 
if other tidal material should be observable (i.e. struc- 
tures associated with extended tidal tails). We have 
not detected such material, to sensitive limits, in our 
Keck/LRIS longslit observations. There are several rea- 
sons, however, why we may expect the non-detection 
of other ex tended emission. Fi rstly, there may be no 
tidal tails. iHibbard et al.l (|200lD present a compilation 
of H I maps of merging galaxies and find 16 examples 
of merging galaxy pairs (of 181 total unusual gas mor- 
phology examples) with an observed bridge and no tidal 
tails. Similarly, the GalM er merger simulation database 
(jChilingarian et ani2010[ ) has many examples of such sys- 
tems. Secondly, the tidal tails may simply not have fallen 
into the longslit of our spectroscopic observations. Fi- 
nally, the tidal material could have too low density to 
render it observable; emission is an n'^ process and a fac- 
tor of 3 less dense gas results ^ 10 times fainter emission. 
We conclude that the non-detection of additional, tidal 
material is consistent with the observed bridge and stan- 
dard merger scenarios. Nevertheless, we are motivated 
to search for additional emission beyond our longslit ob- 
servation, e.g. with an IFU spectrometer. 



Before concluding, we wish to consider a few issues 
related to the 'green valley' nature of the companion 
galaxy. We remind the reader that this is a rare and 
poorly understood phase of galaxy evolution which nev- 
ertheless may be critical to understanding the (expected) 
transition from the blue, star-forming sequence to the 
'red and dead' sample. The first consideration is to un- 
derstand why the galaxy appears in this unusual portion 
of color space. We have found that the peculiar col- 
ors are consistent a declining starburst stellar population 
that would dominate the SED for only another few hun- 
dred Myr. The next obvious consideration is whether the 
galaxy is currently moving toward the blue cloud or red 
sequence. Unfortunately, any such prediction hinges on 
bold extrapolation of the future SFH. If we assume that 
the observed SFR continues indefinitely (i.e. at « 6Mq), 
then we find that the galaxy evolves to becomes slightly 
bluer {u — g decreasing by 0.1 mag) over approximately 
1 Gyr. However if the star formation ceases then the 
stellar populations will passively evolve to become red 
after several hundred Myr. If the galaxies are to un- 
dergo further interactions with each other, then the star 
formation history may follow one of the standard track s 
towards quenching (see discussion in lFaber et al.lf2007l ). 

This system may also provide insights on some of the 
large nu mber of newly disc overed double peaked [O III] 
emitters ( Smith et al.ll2010l) . If observed with sufficiently 
poor spatial resolution, it is possible that the bridge or 
the companion galaxy's emission associated with pho- 
toionization from the quasar may result in the appear- 
ance of a double-peaked [O III] profile separated by ap- 
proximately 20 km/s resulting from just a single active 
SMBH. One could imagine that if caught closer in time 
to a close passage, where the radial velocity difference of 
the two galaxies would be greater, the velocity difference 
may be large enough to produce the ~ 200 km/s offsets 
observed in that sample. This explanation is may be the 
most likely for a subset of their systems (Rosario et al. 
in prep.). 

8. SUMMARY 

In summary, we have compiled physical constraints on 
a galaxy merger at modest redshift which allows for a 
direct test of the merger paradigm developed by state-of- 
the-art simulations. Impressively, the merger character- 
istics are consistent with the standard paradigm for each 
of the comparisons summarized in Table [5l Our anal- 
ysis presents a relatively clear-cut example of a quasar 
being triggered during a first passage interaction. Obser- 
vational characteristics allow it to be directly compared 
with simulations to constrain triggering mechanisms. 

We have copious amounts of information about the 
companion galaxy and obvious questions emerge such as: 
Why is one galaxy in a quasar phase while the other is 
not? Has the companion galaxy already undergone a 
quasar phase? Are the black holes in each galaxy now 
and will they later be on the many black hole-host rela- 
tions? How does a first passage quasar affect the evolu- 
tion of a merger vs. one where neither galaxy goes into 
a quasar phase until final coalescence? Paper II presents 
the insights that we gain from this system about some 
of these questions as well as discusses constraints on the 
quasar's lifetime and isotropy that we get by examining 
the effects of its ionizing radiation. 
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Our work focuses on characterizing the system as well 
as possible. A summary of our derived properties can be 
found in Table |4] and a review of our conclusions follows 
below. 

1. We presented a unique interacting quasar/galaxy 
pair connected by a remarkable 38 hj^ kpc bridge 
of ionized gas visible in emission lines. 

2. While one galaxy is in a quasar phase, its compan- 
ion shows no detectable AGN activity. 

3. After ruling out collisional ionization and radiative 
shocks, we conclude the bridge's source of ioniza- 
tion to be photoionization by the quasar light. 



metallicity is consistent with being a few tenths 
dex smaller than the companion galaxy. 

8. The bridge's velocity dispersion is larger than ex- 
pected 40km/s), possibly suggesting that some 
material in the bridge may be self-gravitating or 
there is some source of turbulence in the bridge. 

9. Using the projected separation, we find the merger 
stage to be between first and second passage. 
This interpretation is consistent with the observed 
bridge and the inferred starburst stellar popula- 
tion. 

10. Based on dynamical measurements and stellar pho- 
tometry, we argue the mass ratio is likely ~1:4. 



4. The companion galaxy's color and spectrum are 
consistent with a declining starburst stellar popula- 
tion. In this scenario, we date the burst stellar pop- 
ulation to be 300 — 800 Myr old. This intermediate- 
age stellar population places the galaxy within the 
'green valley'. 

5. We found the companion galaxy to have a SFR 
of ~ 6 Mq yr~^ and to have a have a gas-phase 
oxygen abundance consistent with solar. 

6. The bridge between the two galaxies was character- 
ized to be dependent on a geometric parameter V, 
which we constrain (see Fig. [5]). We found expres- 
sions for the volume density, column density along 
line of sight, and mass of the bridge as a function 
of this parameter. 

7. With CLOUDY photoionization modeling, we deter- 
mine the ionization parameter and metallicity in 
two separate patches of the bridge. The bridge 
is nearly fully ionized in hydrogen and is primi- 
arly comprised of higher ionization species. The 
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